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The Strains and the Energy in Thin Elastic Shells of 
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Part 1: THEORY 
1. Introduction 


The vast literature existing on the theory of plates and shells can be con- 


4 veniently divided into two groups. The one is concerned with establishing the 


fundamental equations of equilibrium and of motion, the other (which is much 
larger) deals with methods for the solution of these equations for special cases. 


~ We shall briefly review two exponents of the first group to allow later a compa- 


rison with the theory set forth in this paper. 

S. TimosHENKO [1]2) considers an element of the shell. The forces acting 
thereupon are expressed in terms of the deformations, i.e. the extensions and 
changes of curvature, and then subjected to NewTon’s law of motion. The 
changes of curvature are simply defined as the differences of the principal 


1) Presently at Bell Telephone Laboratories, Inc., Murray Hill, N. J. 
2) Numbers in brackets refer to References, page 189. 
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curvatures before and after the deformation. Since the principal directions may 
change by a finite angle, some cases e.g. the spherical shell must be treated 
separately. The effect of the extension of the middle surface on the changes in 
curvature is neglected (page 354). This has strange consequences, some of 
which will be discussed in II.2. It is a disadvantage of the theory that the 
procedure for determining the changes of curvature has to be repeated from its 
beginning for each new case. If the shell is not of a very simple shape, the 
necessary geometrical representation becomes awkward (see e.g. figures 149 or 
176). Moreover, it does not seem very satisfactory that temporarily a shear 
force normal to the middle surface is introduced, on which no boundary con- 
dition may be imposed (page 89). This stems from the fact that the influence 
of this force on the deformation is (and rightly so) neglected. 

A.E.H.LOovE [2] presents an approach which seems to be more analytical; 
however, it is not easy to follow his arguments without drawing figures. In 
fact, the verb ‘to draw’ is repeatedly used. In the general theory (Chapter 
XXIV) nine quantities are used to describe (rather unsymmetrically) metric 
and curvature of the middle surface of the shell, although it is known [3] that 
six quantities suffice for a symmetrical and compact formulation. As coordinate 
lines on the surface the lines of principal curvature are chosen. This special 
choice restricts the applicability of the theory since in many practical cases the 
lines of curvature are not easily found or do not allow to express boundary 
conditions (see II.1 and II.3). In the definition of the changes of curvature 
(Article 324, equations 12, 13) again the influence of the strain in the middle 
surface is neglected. Consequences of this approximation will be discussed in | 
te 2. 

This paper is based entirely on LAGRANGE’S variational principle, which 
provides the most satisfactory approach to both static and dynamic problems. 
The coordinates of configuration and the independent variables may be chosen 
arbitrarily. The dynamic boundary conditions are automatically obtained 
together with the differential equations. If auxiliary kinematic conditions (e.g. 
extension-free deformation) allow no explicit solution for the coordinates, they 
can be handled by the multiplier method of LAGRANGE. 

However, before the techniques of the calculus of variations can be applied, 
the kinetic and potential energies have to be expressed in terms of the coordi- 
nates and independent variables. This is the program for Part I. 

The fundamental idea of the theory of shells is to describe the deformation 
of the entire shell in terms of the deformation of its middle surface alone, thus 
reducing a three-dimensional problem to a two-dimensional one. The necessary 
approximations are these: (1) During the deformation each point of the shell 
remains on its normal, although its distance from the middle surface may 
change. (2) The state of stress is plane in each point of the shell. (3) The strains 
vary linearly across the shell. These approximations are shown to be valid for 


a 
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thin shells, that is for shells the thickness of which is small compared to its 


radii of curvature. More precisely, the errors in the strains are shown to be of 
the order of A/R, where 2h is the thickness of the shell and R the smaller of 
the two principal curvatures. It should be pointed out that the present theory 


] “considers the strains as fundamental and defines the additional curvatures in 
terms of the strains, in contrast to current theories. Approximation (1)—in 
connection with (2)—makes sure that the stress-free plane is parallel to the 
_middle surface. As a consequence of approximations (2) and (3) the state of 


ez 


strain along a normal of the middle surface is described by six quantities. On 


_ the other hand the deformation of a surface in any point is completely deter- 


t 


mined by six quantities. This justifies the approximation from our basic point 


_of view explained above. Corresponding to the strains, there will be stresses 


po 


varying linearly across the shell. They in turn will give rise to three forces 


_(2normal, 1 shear) and to three moments (2 bending, 1 twisting) all parallel 
_to the stress-free plane (Appendix). A shear force normal to this plane cannot 
be accounted for, which is the reason why its introduction leads to difficulties 
when it comes to the formulation of boundary conditions!). External forces 


which may or may not be parallel to the middle surface are represented by 
their work under a general displacement. 


2. Outline of Approach 


The middle surface of the undeformed shell is represented by a vector r which 


is a function of two surface parameters. These define a net of curvilinear coor- 


ra 


dinates on the middle surface. The displacement shall be denoted by the vector 
or. If 9 denotes the density of the shell and 2 its thickness, the kinetic energy”) 


is simply: 
Exn=f [eh (sp or) 4S, (1) 


where the integration has to be extended over the whole middle surface S. 
Both quantities 9 and / may be functions of the position on the surface. The 
expression for the potential energy will be considerably more involved. Assu- 
ming a plane state of stress in each point of the shell, the potential energy [4] 


"per unit volume can be expressed as follows: 


= ay : 
Vie Q py {(Ep» y &qq)" +2 (1 rig LX) (ei ¢ - ene Sea) , (2) 
Y = Younc’s modulus, 
js = Potsson’s ratio, 
e = Strain components (2 &p¢ = Ypq)- 


1) S. THIMOSHENKO, l. c., page 89. 
2) A. E. H. Love, 1. c., page 512. 
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Here the components of the strain, ¢, are taken relative to certain Re | 
p and q, which are normal to each other and parallel to the middle surface. We: 
denote by w the normal distance of a point in the shell from the middle surface. 
The strains, depending on w, shall be approximated by the first two terms of 


their power series in w, that is: i 


|| a {I is 
| En p Enq | | np Ena ; | “pp Xnq | 
ad eh ee a) 


Epa €aa | Enq faa || *pa qq || 


Thus, the e’s and the x’s are defined in terms of the strains, ¢. The é’s represent: 
the strain of the middle surface, the x’s describe the additional curvature due: 
to the deformation. Inserting this into (2) and integrating V across the shell! 
over the range — h to h of the variable, w, we find the potential energy per? 
unit area!): | 
Yh = Belics eee a es . 

eS ye L En + Fad) + otha a) (23 a fgg Sa) } | 

: (4) 

Yh 2 1 apa 4 4 : 
3(0= pL?) { (ny + Aaa) + 2,( — ft) (Wp on a anes Rey eine ; 


The first term, proportional to h, represents, the energy of extension, the 
second part, proportional to h3, is the energy of bending. The total potential 
energy due to deformation then becomes: 


Es -[{ fv oo (5), 


Any system of external loads can be accounted for by expressing their work | 
OW for an arbitrary displacement. Thus, if Z; denote discrete forces at the posi- 
tions r; and 2(w, v) represent distributed forces, )W becomes simply 


OW = )" Z, or; Bie or os. 


Our main task, therefore, is to find suitable directions #, g in every point 
of the shell and express the components of strain along these directions in terms 


of the displacement 6r. 
3. Analytical Representation of a Shell 


Any thin shell is completely described by its thickness, 2 h, and its middle 
surface which shall be represented [3] by the vector 


r(u, v) , (6) 
1) A.E. H. Love, l.c., page 531. 
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which is a function of the two parameters, u,v. While these parameters are 
varying in a region of the (uw, v)-plane, the vector r sweeps over the middle 
“surface (Figure 1). The vectors 


< 


——=r, and ———-=r, (7) 
z ~. 
_are the tangents to the coordinate lines v = constant and u = constant on the 


x 
- 


a 8 ee a 


es \ 'y/ 


Figure 1 


Coordinate system on middle surface. 


surface. The unit vector normal to the surface is 


Ty X Py : (8) 


Now any point within the shell shall be defined by the three coordinates 
(u, v, w) as follows: 


x(u, v,w)=r(u,v,)+wn(u,v), (-ASw =h). (9) 


In calculating strains, distances between infinitely close points must be com- 
pared. Therefore, we need the expression for the linear element. 


tdx|? = (r2+ 2wr,n, + w? nj) du? 
+2(r,r,+ v[r,n, + ron] + wn, N,) du dv 


+ (r2?+2wr,n,+ w?n?) dv? + dw? 


or: : 
|\dx|? = edu®+2fdudv + ¢g dv? + dw?. 
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Hence, the metric tensor underlying (10) has the form: 


ite, ae : 
Lt ¢ 6]. (1) 
lo 0 1| 


A linear element parallel to the middle surface shall be denoted by ds, thus: 


ds = |aX| iy -o - (12). 
Since we consider thin shells, we retain only the linear terms in w, thus: | 
dx? = (r2+2wr,n,) du? | 
+2(r,r,+v[rnn,t+r,n,|)dudv | (13) | 
+ (r2+2wr,n,) dv? + dw?. | 


The omission of the terms containing w? is justified, as will be seen later, if the 
radii of curvature of the middle surface are large compared to the thickness 
of the shell. Using the abbreviations 


nak, 7,.%,=f> r= €. | 


1 (14) 
—r,n,= TL, ees (ny r,n,) = M, rjn,=N, | 
we find 
@=E—2oL, feF—2uWM, g=CGo2om, (15) 
while the normal, n, can be rewritten in the form 
n=(r, x r,] (EG — F2)-3?, (16) 
Here use has been made of the identity 
[a x b]? = a? b? — (ab)?. 
From r,n = r,n = 0 follow by differentiation the identities 
L=-rn,=r,,n, | 
=—T n= rn =r, i. (17) 
=—-rn,=r,,n. | 


The curvature of an arc, r(s), is defined as k — |d*r/ds?|, provided the 
parameter, s, is the length of the arc. Now let the line be traced upon the sur- 
face, r(u, v), and its principal normal, d*r/ds®, coincide with the normal, n, of 
the surface. Its curvature then defines the curvature of the surface in the 


PO EEN eee ee ee 
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direction of the line. Since n (dr/ds) = 0, we find by differentiation: 


d dr dn dr a*r 
ds (n He ee ea al uh) 
and therefore k = | (dn dr)/ds?|. Instead of using the absolute value, the follow- 
ing convention for the sign has been adopted: 


— dn dr 
“aso 


ge (19) 
Thus, & is positive if the normal of the surface and the principal normal of the 
curve point into the same direction. Expressing dn, dr in terms of du, dv we 
obtaia 
dn=n,du+n,dv, (20) 
dr=r,du+r,dv. (21) 


Substituting (20), (21) into (19) and using the abbreviations (14), (17), we find: 


L du2 + 2M dudv-+ N dv? (22) 


a ee es a a 


The maximum and minimum of & represent the principal curvatures of the 
surface. To find the equation for these curvatures we observe that they equal 
the maximum and minimum of the expression 


L du? +2M dudv-+N dv? 


under the condition 
E du2 + 2F du dv + G dv? = constant. 


Therefore, by LAGRANGE’S multiplier method, we have 
(L—AE)du+ (M—AF)dv=0, (23) 
(M —AF)du+(N—AG)dv=0. (24) 


The eigenvectors (du, dv), of this system give the directions of the principal 
curvatures, and the eigenvalues of A determined by the equation 


(TE INS 1G) = (M = AF)? =0 (25) 


are the principal curvatures /;, k,. To prove this, multiply (23) by du, (24) 
by dv, add and solve for A. From (25) follow ViETA’s relations 


EN-—2FM+GL 
hy ks = as a (26) 
LN— M? 
ky Ry be Woxee_ val =e (27) 


160 Ericu S. WEIBEL ZAMP 


If, for a moment, the lines wu = constant and v = constant are chosen as the 
lines of principal curvature, F = 0, M = 0, or using (14) and (17): 


rr, =7,Nn,=T7Tn,= 0. 


This means that r,, and n, as well as r, and n, are parallel. Then the linear 
element (10) becomes: 


| dx? =[]ry|? + 20 |r,| |m,| + w? |n,|?] de? 


+ [| rol? + 2 w|r,| |n,| + w?|n,|?] do? + dw? . 


The terms in w? can be neglected if 


h?|n,|?<2h|r,||n,|, 
that is if 
n<2ttul 2p, 


|r, 
and similarly if 


oR 


where R, and R, are the principal radii of curvature. This justifies the approx- 
imation (13). 


4. Deformation 


Let each point r(w,v) of the middle surface be displaced by a vector 
or(u, v) (Figure 2). Consequently, each normal vector n of this surface will be 
transposed into a new normal n + 6n. Products of deformations will be neglec- 
ted (for instance, from |n+ 6n|=1 follows n6n = 0). We assume that 
during the deformation each point of the shell remains on its normal wits 
distance from the middle surface, however, may change from 


w to [1+e(u,v)]w, 


It will turn out that this lateral strain, ¢, is of no consequence. 
The point, defined by its coordinates u,v, w, therefore, is represented after 
the deformation by (Figure 2) 
x+ 0x=r-+or+(1+e«)w(n+ on). (28) 


The metric tensor, giving the deformed linear element |d(x + 6x)|? will be 
changed accordingly: 


e- de f+ df de. | 
f+oéf gt+dég  d8o¢ 


| O£13 Ofo3 1+ dfs 


b, 
a 


4 


d 
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For the various components we find 
e+ d0e—E>0R—2e(E + OL) (i+), 0g, = @4,, 
f+o0f/=F+ 6F —20(M+6M)(1+e), dg,;=we,, 
g+og=G+0G—2w(N +0N) (l+e), dOf=2¢. 


_ The quantities E + dE, L + OL, etc. are obtained by applying equations (8), 
- (14) and (17) to the vector r+ or instead of r. 


Figure 2 


Displacement. 


To show that ¢ can be neglected everywhere but in 6g33, we use once more 


~ the special coordinates for which F = M = 0. Let us consider the quantity de: 


be =0E —2wdL—2wle. 


Let e, and ¢, denote the strains in the middle surface in the directions wu and v 


respectively. Using a later result, namely e, = 6E/2 E*), we may write de as 


be = 6E (1-w f)—2wol. 


Now 
€ = —m (ey + &) 


and assuming 


[er] = |e 


1) This expression actually equals €,. Since €, will be multiplied by w, the remaining part, 


—x,w, may be neglected as a second order term in w. 


ZAMP VI/11 
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we obtain 

je] S2me,. 
Thus 


é 


& 


<2 and w 2 a Py). 
Af 


Therefore, de takes the form 


de= dE (1-9 4) -2wdL (\9| 2), 
1 


where R, = L/E is the principal curvature in the direction wu. For thin shells, — 


by definition, 4/R, is negligible, that is ¢ is negligible in de. In an analogous 


manner it can be shown that ¢ is negligible in df, dg, 6g,3 and dgs3 so that the 


metric tensor reduces to 
| e+ode f+ 0f 0 | 
|f+ of g+ 6g 0 | 
0 0 L+2¢ 


(29) 


This form finally justifies the assumption of the plane state of stress. The 
components of (29) now are simply 


e+ de=E+6E—2w(L+O6L), | 


{+ 6f=F + 6F —2w(M-+ 6M), (30) 


g+og=G+06dG—2w(N+O6N). | 


This representation of the metric tensor holds also for Jarge deflections 6r as 
long as dE, OF, etc. are small. If they are not, the linear theory of elasticity 
can no longer be applied. 

During the deformation the square of the linear element, ds = [AX eocon ae 
any point in the shell in any direction parallel to the middle surface (du, dv, 0) 
will change to: 


(ds + dds)? = (e + de) du? + 2 (f + df) du dv + (¢ + 6g) dv?. (31) 
From this equation the strain: 
ods 
eater Pe (32) 


in the direction (du, dv, 0) is easily found to be 


E= 


1 bedu® + 2 of du du + dg dv? 
is, 


edu®+2fdudv+gdv2 ° (33) 


Pe eT a 
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By the same reasoning as that used for the curvature, we obtain the equation 
for the principal elongations ¢,, &, 


(de — Zee) (6g — 2g) — (6f —2ef)? =0 (34) 


and by VietTA’s relations: 
_ €0g— 2f of + g de 
sig, Sie Fis gen pe)" 
be dg — Of? 
id UN se is 


5. The Components of Strain 


The expression (33) gives the strain at any point (w, v, w) in any given 
direction (du, dv) parallel to the stressfree middle surface. We shall choose in 


~ each point of the shell orthogonal axes, # and q, parallel to the middle surface. 


With respect to these axes we shall express the components of the strain 
tensor: 
(37) 


i Epp €pa 
| 
I} 


Epa faa 


This reduced two-dimensional form is sufficient, since it completely determines 
the plane state of stress. Once the directions of the p- and g-axes are given in 
terms of increments (du, dv) the extensions ¢,, and ¢,, are readily expressed as 
the strains along the axes # and gq. To find the distortions é,,, the strains are 
calculated along the axes of a new coordinate system %, (Figure 3) which is 


Figure 3 
Coordinate transformation in a point P. 
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rotated by — 2/4 with respect to the (, ¢g)-system, so that: 


p 1 1//]| +] 
ae Hi |. (38) 
| V2 | | 
ca a ALIN 
The old components in terms of the new ones are 
| to ft 1 
ll En» epg | |= (Snot pa) oe ong — > (Exe — Evy) 
| Sit : be (39) 
| | Vv iv | 
1€pq €aq. | aes (ei acaas @ (Exe t Evy) — Exy | 
1 | | 
in particular 
1 
Spa. = iB (Eyy > Ean) - (40) 


To obtain symmetrical formulas we choose as x- and y-axes the bisectors 
of the curvilinear axes, u, v. This defines also the #- and g-axes. The transfor- 
mation between the w-, v- and the x-, y-axes is!): 


awe (2V/esin-$ By" ; (2V/e cos sy) | | de || 
| = AS Hil (41) 
de 1 — (2Vgsin$ ay : (2Vgcos2) | dy 
where 
cosw = f(e g)~*”. (42) 


It is understood that the positive values of the roots are taken. The linear 
element simply becomes: 


ds? = dx? + dy?, (43) 
The transformation (41) is a function of u, v, and w and is in general not in- 


tegrable. This transformation makes it possible to express the directions of 
the x- and y-axes in the variables u, v as follows: 


x-axis: { du, dvs, = { gn. ee be (44) 
y-axis: {du,dv},u{e?, p-Vet (45) 


The symbol ~ stands for ‘is proportional to’: we are only interested in the 


1) This rranctorawtien may be found analytically, or directly from Figure 3. 


ee AY A ee 


ee ee ee ee ee 


c 
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the inverse transformation of (38), namely: 
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directions of the axes, not in the length of the vectors representing them. Using 


together with (41), we find also the directions of the p- and q-axes: 


a Bi ag me En/2 a ‘ QO na SP @ takes wo 
p-axis: {du, dv}, © le (cos > + sin ar g (cos pha eie ale 


g-axis: {du, dv}, » je" (cos. — sin), goes (cos + sin 5). 


(46) 


(47) 


The p- and qg-directions can be brought into a more convenient form by using 
equation (42): 


{ du, dv}, Aine 


oa 
_ 
3 
—i 
= 
+ 
—— 
aay 
| 


ale are [1 = (1 4 ere 


= 1/2) . (49) 


By comparison of (48), (49) with (46), (47) for small values of f/(e g) 1t follows 


that [1+ 


(1 — f2/e g)/7]"? is always positive, whereas [1 — ( 


1 — fre gy?” 


takes the sign of /. Upon substitution of the values of du, dv for the various 


sions for the strains in these directions are obtained: 


Eq 


pees Bice) ote evel ih ei eg 


oe 4 (eg — f) 
[g de + 2 (eg ne df +e dg] [1 — f (e g)11 
4 (eg — f?) 
ee seg tt Peg 
4(e¢ — f*) 


Loe fh (0 — fre g)) — 2 f of +e dg (1 + Gye | 


rice = -) 


Using equation (40) we readily obtain the shear strain €y, 


. f (gle)? eae of + f (e/g)! dg 
ae 4 (eg — 7?) 


directions x, y, f, 7, into the formula (33) for the strain, the following expres- 


(50) 


(51) 


(54) 
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ViETa’s relations (35), (36) provide a welcome means to check the strains, 
since they give the two invariants of the strain tensor: 


€1 4 Eo = Epp + Egg = Exa tT Evy » 


pre 2 Lo 
Ey En = Epp Eqgq — Epa = Ex Eyy — Exy> 
where 


ut 
oO Mane (xe = fo) : 


The components of the strain tensor as given by (52), (53), and (54) are not 
yet in the required form. According to the program outlined in the introduction, 
these strains must be linearized as indicated by equation (3). From (15) and 
(30) we obtain: 


d¢e=OE—2wdL, df=dF —2wdM, dg=dG—2wdN. (55) 
Substituting (15) and (55) for e, f, g, de, 6f, dg into the expressions for the 
strains, we arrive at functions in w; for instance (54) becomes: 


5 a (E—2wL) (G —2w N) —-(F—2wM)?]"* 


x {(F —2w M) (G — 2w N)!? (E — 2w L)-?” (6E — 2 w OL) 
— 2(E —2wL)” (G — 2w N)? (6F — 2w 6M) 
+ (F —2wM) (E—2wL)'" (G—2wN)-'? (6G —2w6N)}. 


Similar expressions are obtained for ¢, ,, &,,- From these expressions, the coeff.- 
cients of the first two terms of their power series 


e=E—xW (3) 


have to be extracted. Since this process is rather lengthy and offers no particu- 
lar interest, the results alone are given: 

The strains of the middle surface are: 
5 — GOE (1+ [1— FYE G}") — 2 FOF + E£6G (1 —[1— FYE G}) 


pp 4 (EG = yan Se ae aera , (56. 1) 


- _—s«-F (G/E)* 6E — 2(E sal 2 OF + F (B/G) 6G 
Ena eg 3 Me. a7 (E =F) os a , (56. 2) 


2) Gen t=] ae Oe ee 1+ [1— FYE Gp" 
€qq 4(E Cs F?) [ | iE Z . (56.3) 


j 
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- 
‘The additional curvatures!) are: 


2 — 2 (ky + Re) Epp + Pes 

P «(2-9 0-Ay "re G-28 +8) | 

+2(N6E+G6L) [1+ bes | eer | 57.1) 

i +2 (L6G + E 6N) [1 - (1-+<) |}: | 

Be, = — 2+) 

ee Pe ye | 

pag. ae pe Pere) 

+t a(b ag) | 

=| 57.3) 
‘| 


E 
2 (L8G +E 6N) [1 eee) 7 4 (M 6F + F 6M) 
ee 


The terms proportional to (k, + k2)1) in the equations (57) arise from the 
linearization of the denominator (e g — f?): 


aa eo 1 EN-FM+GL 
FE a a oe |, 


which by (26) reduces to: 
(eg —f)-1 = (EG=— F2)-1[1 + 2 (hy + fy). 


The equations (56) to (57) give the desired strains and additional curva- 
tures due to deformation as functions of the displacement 6r, since all quanti- 
ties 6E, dL, etc. are functions of dr. The formulas are simplified considerably 


1) See last paragraph of this section, page 168. 
2) These are the terms which are neglected by A. E. H. Love, S. TIMOSHENKO; see introduction. 
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if the system of coordinates on the surface is rectangular; then F = 0. If in 
addition, the coordinate lines are chosen to be the lines of principal curvatures, 
M will disappear. Two important special cases are listed below. 

(a) Lines of Principal Curvatures as Coordinate Lines. F = M = 0 


4 jet ts or et = 6G 


eae In ayeG" = 26 ee 
OL LOE ON N 6G 
lr ee aS ee ae Pee, 
(57a) 
a es e ‘N), OF | 
Xnq — WEt G = a0 \e 7) VEG : 


(0) Average Curvature k, + ky = 0, Orthogonal lines of Zero Curvature as 
coordinate Lines. F = 0, L=N=0O0 


ie A ae me) it os 


Epp 2E? Ena 2VEG’ Eqq 26 ; 
OL VI ON M 
ino Ee EEO ee a we oe | 
(57b) 
da Mat Page ee , OE | 
Vee 2 (ee eae 


It may be noted that in the first case x, ,, x,q and in the second case x, , can be 
expressed as total differentials: 


N 


ik 
(4) Xp» =O, = Ohy, Haq = OG 


=e 


M 
(0) %ap = VEG. 


Only in these instances the x’s merit their name ‘additional curvature’. 


Appendix 


Although we do not make use of the stresses, stress-resultants and stress: 
couples, these quantities will be given here. According to the wellknown rela- 
tions between stress and strain we find from equation (3): 


re 


hs z ba 
Ope Tis pe (Epp + [ Eqq) — W 1 2 (pp + Mt Hea) » 
— ‘4 ra e s | 
or anal 7wr we (Eqq + Epp) — @ Tage (“aa + M Xp») » 
Gino pad 
2a = A Be as pa Ww 1 cm u Xnq 


(6b) 


ul wie) 
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By definition!) the stress-resultants and the stress-couples are given by 


+h +h 
Ky, =| on» dw, Myy = —Mga = — | Op 0 dw, 
—h —h 
+h +h 
Kea = | O40 de Mog = — | O00 dw , 
—h -h 
+h ih 
Kya = Kyo = | O50 dw; Me, = Op, w dw. 
Sh 27 


- Carrying out the integrations one obtains: 


4 


The stress-resultants 


Kop = pat (Epp + M €gq) (normal force), 
1 Ca an (aq + Epp) (normal force), 
i = on de (shear force). 
The stress-couples 
Mo5= - aa 5 og = —M aa (twisting ee ctgent: 
Mes 3 a “=i (%aq + { %pp) (bending moment), 
My = — 5 ai (pp +  %qq) (bending moment). 


Thus the stresses, stress-resultants and stress-couples are expressed in terms 
of the strains. of the middle surface and of the additional curvatures. The 
names customarily used for these quantities are indicated in the parenthesis. 


Part II: APPLICATIONS 


1. Plate in Oblique Rectilinear Coordinates 


Although this simple case can easily by treated in a conventional manner [5] 
it will serve here as an illustration of the method. 

Let the vectors 1, jf, k denote the basis of a cartesian coordinate system 
(Figure 4). The (i, j)-plane which shall be the middle surface of the plate can 


1) A. E. H. Love, 1.c., page 461. 
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be represented parametrically by the vector 


r=(u+vcosa)i+vsinas, 
or symbolically 
r={u+dcosa, vsina, 0}. 


Here the parameters u and v are oblique coordinates in the middle plane of the 
shell. 


Figure 4 


Oblique coordinates for plate. 


If only deflections, ¢, normal to the plate are considered, the displacement 
field takes the form: 
dr ={ 0, 0, C(«, o)}. 


so that for the deformed middle surface this representation holds: 
r+ or ={u + v cosa, vsina, C(u, v) } : 


The functions E + dE, ... are calculated by applying the formulas (I. 8, 14) 
on r + or. The tangent vectors are: 


(r+ or), ={1,0, Cu}, (r+ dr), ={cosa, sina, ¢,}. 


Neglecting quadratic terms in the displacement we find: 


E+ 6E = 
Since 
6b =0, for Cs-0, 
we have 
pg Wt lg 
Similarly 


a a 


la 


a AY) yea 


Vol. VI, 1955 The Strains and the Energy in Thin Elastic Shells yal 


The formulas (I. 56) show that €, , = €) 4 = &yq = 0. The deformation, therefore, 


is free of extension. The normal of the middle surface is parallel to the cross- 
product of the tangent vectors: 


n+ 6n=A(r + or), X (r+ on), 
=Asina{—¢,, — (sina)-1¢, + ctga ¢,, is 


The multiplier 4 is determined by the condition |n + 6n| = 1: To the first 
order in the displacement, A sin « must equal 1. 
Next we find the second partial derivatives of the vector r+ or: 


(r+ Or), ={ 0, 0, ee) 

(=) OF) ..6 Rat aes 

ty Or), = 10, 0, foe } - 
Formulas (I. 17) yield L + 6L = (n+ On) (r+ Or)uu = Suu 


Since dL = 0 for ¢ = 0, we infer 
ba), 0b = Cu 


and similarly 
M=0, 0M=,,; N=0, ON-= Cie; 


Applying the formulas (I. 57) we obtain the curvatures: 


1+ sing . 1 — sine c COS a c 
>» = So ---—— — 
BD Sie ee 2sin2?a °°” Sthaletose eee? 
1 cosa rs 
%nq = “and = og a an Raia ar em 
* sin*« 2 sin*a 
1+ sing t 1 =—sing . GOS) 
%, = = #8) 0 
ae 2sinton "°° " Zsin®a °** sin2?« >“? 


Substituting these expressions into equation (I, 4) we arrive at the potential 
energy per unit surface of the plate: 

YR? 
3(1-—u 2) \sinta Tate 
This energy density is identical with the one found by H. FAvRrE [5] by means 
of a coordinate transformation. 


2(1— pH) ree ae I. 


oz [can't Cog — 4 COSe Curl? + sin2a buv Suu svv 


2. Cylindrical Shell 


The middle surface of a cylindrical shell of radius, a, can be represented 


parametrically as follows: 


r—acospi+asingj +zk={acosg, asing, ae 
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where @ and z are the surface parameters. Let the displacement of a point of | 


the middle surface be described by the vector 
or ={é cosy — n sing, & sing + y cos, C}, 


where §, 7, €¢ are the normal, azimuthal and longitudinal components of the 
displacement. The vector 


r+ or ={ (a+ ) cosp — nsing, (a + &) sing + n cosy, z+ C} 


represents the deformed middle surface. Putting a + & = 7, the first derivatives 
become: 


(r+ or), a (7, — 7) cosy — (r + y,) sing, (r + ,) cosp + (7, — 9) sing, ay 
(r + or), ={7, cosy — 7, sing, 7, sing + , cosy, 1+ ¢,}. 
The second derivatives are: 
(r+ Or), ={ (Yop —7 — 2%) COSP — (tp —N + 27,) sing, 
t.o—? ~2y,)sing=- OL, = 942 7) COsg, Lays 
(r si One == { tae > Nz) COS ong (7, = Ny 2) sing, (r, “TF Hiss) COSY 
+h ig ay) sin piggy t9 
(r+ or),,= or COSY — zz SING, 7,, SING + H,, COS G, Get 


The various scalar products which now have to be formed are evidently 
invariant under a change of the origin of , that is under the transformation 


P=P+ HP, Or'(y’, 2) = Or(g, z). 


Therefore, we may put m= 0 and thus considerably simplify the algebraic 
operations. 


We find according to equations (I. 14) 
(r+ or)i=E+ dE = (€,—m)?+ (a+ét Holme oe 
Retaining only the terms which are linear in the displacement, we arrive at: 


E+ dE =@+2a&+2an,. 
Similarly . 


E+ oF =an te G, Gt OG 1 ei 2 Ga. 
Hence 


B= @, OF = 2a(€+ ise PaO}, aan ee eCele oC ore 


AVR ee eae Oe eee 
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The normal becomes after aes and linearizing: 
non = {i es , 8}. 
Applying the formulas (I. 17) we find: 
L+é6L=——-a—£+7,,—2y,, M+6M=&,,—,,. N+ ON =,. 


Hence 


L =—a, 6L =§,,— 21, —§ 
MM = 0", 6M = €,,— Ns 
Sela ON Se: 


Now the formulas (I. 56a) can be used to obtain the strains 


& 
S ae a np = Nz lop S = 
fog 4 —) Epz me ae tore 


and the formulas (I. 57a) furnish the additional curvatures: 


y & ets A Eg Zz Nz Ce 4 a 
: toz = =. M4 é 
sage a? Dg a 2G Dir, Wake BEC 


For displacements which produce no strain in the middle surface, 


Lng = ton = as = 0 
requires 
= Gn et 
f=—n,, m=——2%, &=0 
and therefore 
N22 oe = 0. 

Thus the additional curvatures become: 

£ hy ain Pes 

UE ies cae ees rape Ot 


PP a a 


These formulas are identical with those found in Literature (e.g. A.E.H.LovE, 


ie, Article 319). 

In the general case, however, t 
since previous authors neglected t 
face on the curvature. A.E.H.Love gives 

curvature the expressions (in our notation) : 


here is a discrepancy which is to be expected 
he influence of the strain in the middle sur- 
(l.c., Article 334) for the additional 
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If we form the differences x — x’, we find indeed, 


yi St Ne. eve G "ez 
*op — *pp a, a ee 


Hzz Hag =) 
that is, the differences are proportional to the strains of the middle surface. 


(Since k, = 0, x,, = %;,). We can illustrate very clearly the effect of neglecting 
the strains in two simple examples. 


(1) Let €=constant, 7=€=0: 
Our equations give 


= ae 2 Z ‘3 
[ag = Bo Egz = &22=0, gigi 


It is quite clear that a dilation, €, of a cylinder changes its radius and therefore 
its curvature, so that %p, Must not vanish. 


(2) Let - n=bsing, €=C=0. 
Our equations give 

= b = 

Eng = 7 COSY, Sr =—&e=0, w= % 


Standard textbooks give the same except for Hogs 


, b 
%og = sr COSY . 
It is evident that there cannot be a change of curvature associated with the 
assumed displacement and x,,, should be zero. 

It is true, however, that the terms discussed here are usually insignificant 
in their contribution to the potential energy. 


3. Helical Shell 


The middle surface of this shell is generated by a straight line, normal to 
a fixed axis, which rotates about the axis as it progresses along the axis, the 
angular rotation being proportional to the axial displacement (Figure 5). The 
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parametric representation of the surface is: 


Years Bea ae | me 
r=xcos--t+xsn7- jt yR= {xcos 2, “sin, vy}. (1) 
x and y being the parameters. The lines y = constant are straight, representing 
the generatrix in different positions. The lines x = constant are coaxial helices. 
_ The distance 22a, known as the pitch, represents the displacement of the 
generatrix for one revolution. For instance, the middle surface of a propeller 


blade may be a section of such a helical surface (Figure 7). 


; 


Figure 5 


Coordinate systems used for helical shell. 


The general displacement or of the shell shall be defined by the three 
functions &(x,¥), (x, ¥), ¢(*,¥): 


ér=e,f£+e7+e,¢, (2) 


where @,, €,, €; stand for the unit orthogonal vectors 


e,= | cos, sin=., ae | 
eu= — sin - 7 C05 x, of, | (3) 
é,==1.0,. 0; 1}=k. 
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In terms of these unit vectors, the deformed surface is given by . 
r+ or=(x+é)e, + e.+ (y+ ¢)e;- (4) 


The tangent vectors are (parameters as indices denote differentiation) : 


(r+ or), = (Ls Sa) €, + Hz @p + Cn 5, (5) 
(r+ or)y= (E— let (m+Z+slet+(lt+iye- ©) 
Here use has been made of the identities 
0G. eee de, 1 
Ge tn ee eg ere (7) 
From equation (5), (6) using (I, 14) one obtains: 
E=1, F=0, G=1+ (2). (8) 
and 
Oe iF ee 
il 
OP =| an ae 2b, 4 aC), (9). 


. 


2 fa 


iG 2 [els f) bat) 


In general the bending energy in thin shells, being proportional to h3, is much 
smaller than the extensional energy which is proportional to #. Thus the fun- 
damental mode of vibration will, in general, be a pure bending mode!). We 
shall therefore first investigate inextensional displacements by imposing the 
conditions dE = 6F = 6G = 0: 


Lo — OE ae, + at, =0, (11) 
x § ‘ 
= (ny + =) + Cy=0. (12) 


The general solution of this linear system may be obtained in the following 
way: Differentiation of (11) with respect to x, according to (10), leads to 


Rig re Coe ee) (13) 


Since singularities on the axis (x = 0), cannot be tolerated, the functions EtG 


1) Lord Rayreicu, 1. c., Article 235b. 


q 
4 


. 
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_ can be represented by their power series in x 


4 ‘ 


~ 
i 


co 


n= +h Vay) 2, (14) 
c= 7 Wal) 2 (15) 


Substituting these series into equation (13) yields: 


co co 
>, Vania — 1) ge Saw, m (a= 1) pains ca Ue 


_ Therefore, V5, V4, Wo, W,, remain undetermined, but W, = 0 and 


1 w= 1 
va, Vaya? fon 222 (16) 
Thus the two functions take the forms: 
n=At+Bu+ D'V,ly) x, (17) 
p=C+Dx—-% J) Vaal) ee, (18) 


pier ALB. C/D, V,, (2 2)_are still undetermined functions of y. After 
substitution of these expressions into equations (11) and (12), we solve them 


for &, and é: 
t= 2—D, (19) 


a 


-m 


a 2 
ibaa Rar C,+ 4 (A, + 4D) +aBx+ 3 av ny apa * 


Since, according to (10), € is independent of x, C,, B, and (V,)y (7 2 2) must 
be zero, that is: C, B, V, (nw = 2) are constants and 


ge Ds (20) 


J 
a? a 


Differentiating (19) with respect to y, we find: 


A 
Re ere ae (21) 
From the last two equations one derives: 
pe = 2D, (22) 
s 2 
a “w= a A, . (23) 


ZAMP VI/12 
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Let ®(y) be an arbitrary function and put 


2 , \a 
aac (24) 


the prime denoting differentiation d/dy. Upon integration of (22) and (23) we | 


obtain: 
A i OF oats, (25) 
fe ie eae | (26) 


The substitution of (25) and (26) into equation (19) shows that the integration 
constants a, 6 must be equal, so that they may be absorbed into the function @. 
This disposes of that part of the solution which is linear in x. Now consider the 
part involving the series. To satisfy (16), we put: 


Va = —Anyilm+ 1), (27) 
W,, = a-1 A, (nm — 2), (28) 


where the A, are the coefficients of the power series of some regular function 
le, ¢) 
F(a) = pr ARIS s 
0 


It is easily seen that the V,, are the negative coefficients of the series a dF |dx 
while the W,, are obtained as coefficients of x3 d(F(x)/x*)/dx. Therefore, the 
following representation holds: 


ood yin af(x : 
25 oan a mcs (29) 


where 


L eo gen 
ye FN) 


in order to obtain the correct lower summation limits. (The prime is here and 
subsequently used to denote differentiation of any function of one variable.) 
Thus f(x) may be any regular function for which 


P"(0) = #0) = 70) =0 . (31) 


: 
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ye eet ee 


_ Assembling the quantities given by (25), (26), (29), and (30), the solution (17), 
_ (18) appears in the final form: 


3 f= 9, (32) 
} n=—(0"-“)+Bx-f, (33) 
¢== (o" + 5) oe (4). (34) 


_where ®(y) and f(x) are arbitrary functions. 
If f does not satisfy the conditions (31) but only the one that f(x) be regular 
for x = 0, the equations (32) to (34) still give a valid solution. However, the 
terms in { may then also contain the linear functions provided already by the 
other parts of the solution, so that all the terms are not necessarily independent. 
For instance, choosing 


fay act At n+ 5 Be “PS PSCH0. 
where A*, B*, C* are constants, the solution becomes: 
a), Ho A* + Bex, C2C*4 At. 


_ which is the same as if one took: 
S2=Atg® B=Bt, C= C*, f= 0. 


The solution (32) to (34) represents all displacements which result in pure 
bending. They must contain the trivial cases of rigid translation and rigid 
rotation. Let us check some special cases: 

Translation along the k-axis: 


B=@=/f=0, C+ 0. 
Rotation about the R-axis: 
C=@G=f/=0, B+ 0. 
Translation normal to the R-axis: 
B=C=—j=0, @ satisfying ©" + ue =. (). 
Rotation about the i-axis: 


; Tigray y 
Bato pe OG C= —- Sint COS 


a 


7] 
| 
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We proceed to calculate the additional curvatures associated with the 
displacements (32) to (34). Since rigid displacements which produce no strain, _ 
offer no interest, we put B= C=0. The tangent vectors and the normal — 
vector for the deformed surface must now be found. The tangent vectors (5) _ 
and (6) become: . 


(7+ or). =e, — fre, + [2 (0+ 2) 427-17] e,,~ Gf 
(r+ dr) = [—7 (@"+ 2) += p]e+[- (0+ 2) +2]Je, 


+[1+ 2 (or+ 2)]e ie 
a ae oa 
The normal is given by equation (I, 16): 
n+ 6n=[(r+ or), x (r+ 6r),|[(E + 6E) (G+ 6G) — (F + 6F)?] 17. (37) | 
The displacements considered are free of extension; therefore, 
OE = 6F = 6G =0. 
Retaining only linear terms in the functions ®, / we have: 


nsdna([-(1+2)rar- a (0+ S)] 0 


a> a 


+[a-Ee Bye 09) 


+[2-(0"+ 8) Jed (a+ y™. 


The second derivatives of the vector r + ér appear as: 


(r+ 6r)ng=—f" ey + — fey, (39) 

(r+ Orjsy = fe +7 et — (6 + oz) &, (40) 
(P+ bry = — Ze ert [ae — (Sa +z) B) es | 

+ = (6" + <)' e, | 


Now equations (1,17) can be used to give the functions L + 6L,M+6M,N+6N. 
All terms of second order in @(y) and /(x) are neglected. The terms containing @ 
or f as a factor represent the increments OL, 6M, ON, so that L, 6L, M, 6M, 


| 


; 
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_N, 6N can be found individually: 


L=0, M=—(a2+%)"'?, N=0, (42) 


2 


éL = (1+ ew 


a2 


ia fj”, 6M=0, 


év=(1+ 44)" (Se +z) +2, (x f"—f)]. | 


a2 a 


_ Equations (42) express that the parameter lines are not lines of principal curva- 


ture, but rather the lines of zero curvature. From equations (I. 570) the addi- 


- tional curvatures are obtained. 


oe (1 i za ra ren i 
ete 5 (44) 
Ayy — (1 + =) a [ee is ar) Pa as (x f" = /)| - 


a 


If a increases indefinitely, the helical surface becomes a plane. Making the 
substitutions: 


i > fe! > D,, (45) 
and letting a > oo, the additional curvatures become: 
apt, Har =0, My =O" (46) 
while the displacement is simplified to: 
£=0, = —ix)- Oy), =0. (47) 


This displacement is normal to the plane, as expected, and for the curva- 
tures the familiar second derivatives of the displacement are obtained. Equa- 
tions (47) do not represent the most general extension-free displacement of a 
plate, since the strong restraints imposed by the helical shape of the shell do 
not gradually relax; they persist, regardless how large a is made. 

This limiting process suggests that for a long narrow shell extending along 
the R-axis (Figure 6) the function ®(y) will suffice to express its bending. On 
the other hand, the function f(x) will describe the bending of a shell which 
extends far along a generatrix (y = constant) of the helical surface such as a 
propeller blade (Figure 7). 

‘According to (I.4) the energy of bending per unit surface in the first case 
(shell extended mainly along axis, Figure 6) is: 


Vosges wer lat a) 9 ) 
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Now let the function ® be not only dependent on v but also on the time, ¢, and 
denote the derivatives by: 


Figure 6 


Helical middle surface. 


Then the kinetic energy per unit surface according to (I.1) and (32) to (34) 
becomes: 


2 


T a E (0')2 + a2 (o" og i + x2 (0 + aa : (49) 


In the case of a propeller blade, extending mainly along a generatrix, 
(Figure 7) we describe the bending by the function f(x) only: 


EQ. et ae C= +4, 


Here the potential energy per unit surface becomes: 


= yon 1 1 Ks t ast 
CS aie leit ae Lae taste tect 


—2(1~ wf" (eff), 
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and the kinetic energy appears as: 


T=oh|(/)?+ (af — Ne 


a2 


Figure 7 
Propeller blade with helical middle surface. 


Thus the potential and kinetic energies are given for both cases in terms of one 
arbitrary function which describes the displacement. Together with suitable 
kinematic boundary conditions these energy-expressions define completely 
the mechanical system. 

Let us briefly discuss a numerical example. Consider the helical shell of 
constant width (Figure 6), which is built in at both straight edges (y = + /) 
but free at the helical edges (x = + 0). We wish to find the fundamental 
natural frequency as a function of the angle, 2 m, between top and bottom 
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edge. Clearly 1 
ey 0 i} 
are (50) | 
The end conditions are: | 
or =0 | 
F) fory=+/. (51). 
Oy (n or) =, 
Therefore we have 
f=n= l= 0, for y=+1, (52) 
and 
On eo ee : e 
=; ata ee for all x and y= +1 
or 
dn Of 


Let us first assume pure bending as described by equations (32) to (34) with 
B=C=f=0. Equations (52) and (53) then require 


P= P= OG" =O" 20, for “y=+1. (54) 


Since there are already eight kinematic boundary conditions there will be no 
dynamic ones. As a4 ->oo.and y +0 the shell becomes a flat plate and its 
deflection is essentially given by 7 = ®”. From the conditions (54) follows, by 
applying twice the mean value theorem, that ®” has at least two zeros between 
— 1 and /'). The fundamental mode of the shell will therefore not converge to 
the fundamental mode of the plate (nor to any other as can be shown). 

This discontinuity stems from the condition of pure bending é,,, = &,4¢= gq = 0- 
For a deflection corresponding to the fundamental mode of the plate, the middle 
surface of the helical shell experiences shear strains, ¢,,, proportional to the 
angle @, as will be shown below. Such strains, however, have been excluded 
through the conditions (10), (11), and (12) of inextension. 

The extensional energy associated with this strain disappears continuously 
as m tends to zero, thus rendering the argument for pure bending invalid. 
Physically speaking the condition of inextension is not an exact one as given 
by the equations (10), (11), and (12), but is approximately imposed by the 
generally high energy associated with extension?). As @ tends to zero this 
physical constraint gradually relaxes unlike the mathematical conditions (10), 
(11), and (12). 

To represent this situation we have to give the system more freedom by 
introducing a new coordinate of configuration. The bending of the flat plate 


1) This FO not be true, if the shell were clamped at only,one straight edge: (54) would 
apply only to y = —J/ and ©" could be made to have no zeros between —/ and I. 
*) Lord Ray.rian, l.c., Article 235b. 


; 
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_ may be described by 7 = — H(y), § = € = 0. If we make the assumption that 
for a slightly twisted plate the displacement is still normal to its middle 
- surface, we are led to 


x 


f= 0), Wastes Ga els 


Superimposing this displacement onto the inextensional one we arrive at 


g=-29', n=-(9"-= +H), c=*(0"'+4+H). 65) 


To make @(y) and H(y) independent, we require that they satisfy indivi- 
- dually the boundary conditions (51). Thus we find 


Vem S00. for »si, (56) 


_as boundary condition for H, (54) remaining valid. 
As will be seen later, H will describe the deflection for small angles, whereas 
@ takes over for large ones. In the transition both contribute significantly to 
_ the deflection. 
. The next step is to derive in much the same way as before the quantities 
_6E, oF, 6G, and 6L, 6M, dN. Then applying equations (I. 56b) and (I. 57b) 
_ we find the strains: 


big ete (1+ ay" 2 ig =0 
and additional see 
pS pee we (1 4+ a se 
win= eB)" [ide + aos eB), 


The terms containing H/a? in the additional curvatures will be neglected, since 
- the energy arising from them is about 9?” h?/l? times smaller than the energy 
arising from shear. Substituting these values into equations (I.4) and (I. 5) 


we find 
+1 +b 
im Ae ae 1 aes 
$2 =| Hi +p 1+ (x/a)? a? 
=1 -b 
Y hs 1 emacs las ae 
eee : = _ He ds 
aaa Te (ie? Lage ae) + | as, 
where 


as =VEG—F* dvdy=|/1+ (=) dx dy . 
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Making the substitutions 


Halk, Dale. yoo p, 


the potential energy can be given the form 


Y h3b pb 
Eyot = 3 (1 — p?) 1 A ( 1 ) 


[loam Exes (ee of weed a. 


-1 


+1 (57) | 


where 
x 


2 da | 
A(x) = ahrewnice (58) 
0 


Indicating the time derivative by a dot, we obtain from (I.1) and (55) the © 
kinetic energy 


Exin = QhOE t (72) [tee =P P+ K)? + 4 g? (22) ap 


; (59) 
oe oe : ; . 
+ (2°) | ("+ + Kap |, 
where e 
B(x) = 2 7 (1+ a2)! do , (60) 
Gla = zie (1 + a2)" de. (61) 


The energy functions (57) and (59) allow the application of the Rayleigh- 
Ritz method to compute approximately the fundamental frequency, », for 
various angles y. To make the problem definite the following ratios have been 
chosen: 


is 
aes A, pe ole, w=03. 


For the functions Y and K the simplest polynominals, satisfying the boundary 
conditions (54) and (56) were taken 


PF =A, (6? —1)*, K=A, (62-1). 


Br en 


=~. 
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Thus the variational problem is reduced to a minimum problem in two variables 
Ay, A, (with a auxiliary condition). 

_ The accuracy of the frequency depends on how good a guess is made on 
the shape of the functions Y and K. Using more than two indeterminate 
variables would improve the results, but the numerical labor involved already 
in the two-term approximation discourages any such attempt. Automatic 
‘computation is the only answer to the demand for greater accuracy. 

The behavior of the solution shown in Figure 8 is easily understood by 

examining the characteristic equation and its dependence on ¢. For sinusoidal 
motion of the frequency, v, the maximum energies will be of the form 


E pot max = sa aT [Yas Ay + 2 Vig Ay Ay + (Vs) dey az V8) ae | (62) 


Exin max = ohbis (2% v)?(T,, At 4 Ty. A, Ay + The a5]. 


eee Pea ee 
wes \ batho! 
25.0 i A 
oy Matches 
20.0 \ 
SS 
> 17.5 ‘ 
6 \ 
Zz 15.0 
Ww 
= 
Ww 12.5 
- [toed sewed pce aed 
7.5} 
9 
‘ a ee 
“ Pe eet 
0 
0 AoeeGGn eG «120.«(«1S0~”*«BO.:*:*=C RTO 240270 


ANGLE, , IN DEGREES 


. Figure 8 
Resonance frequency v, versus twisting angle, 29. 
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with V;, and T;;, depending on g. Using the abbreviation 


es! |/ i h 
MO PS ayo te 


the characteristic equation may be written as follows: 


Va — (p-) Ta Vis (EY) Tia 


Vin —($-) Ta. VO + 0% 55 VO - (2) Te 


Vo 


Its solution is represented as a solid line in Figure 8, which shows the ratio. 
»/¥ as a function of y. (More points of the frequency curve beyond 180° have | 
been computed indicating an increase in frequency for larger y values; however, | 
it is doubtful whether the function Y for these values of @ still gives a fair 
approximation to the true deflection, so that the v/v) values may become in- : 


creasingly inaccurate.) 

To obtain a clear understanding of Figure 8 we observe that equation (63) 
has the form of the characteristic equation of two coupled oscillators, whose 
masses, spring and coupling constants are functions of g. Oscillator 1 represents 
the extension free bending generated by YW, while oscillator 2 represents the 
deflection normal to the shell described by K, which gives rise to bending 
energy V{)) and extensional energy g? (/2/h2) be 


Let us, for a moment, consider the oscillators decoupled and let »y and vg 


be their individual frequencies. If either of the two becomes very much larger 
than the other, the fundamental frequency of the coupled system will 
approximate the lower of the two. 

If the shell were infinitely stiff against extension (V3 = co) its fundamental! 
frequency would be given by 


a 
m= Vas 
which is represented as a dotted line in Figure 8. 

For small angles, »< hil, the spring constant Vj, is much larger than 
Va + gy? (12/h2) V&, since the elastic line has two nodes for inextensional 


bending (page 53) but none for the deflection K. Therefore the shell will 
vibrate in the mode given by K at the frequency 


As tends to zero the continuous approach to the correct frequency is ensured 
since K describes the fundamental mode of the flat plate. For increasing y the 
stiffness against extension V2) y2 12/n2 grows rapidly (/?/h? = 130,000) and 


— 


ee 
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4 


— 


. 


without limit. This causes yx to rise sharply. As vg becomes comparable to 
Vy, both modes, ¥Y and K are simultaneously excited and vy passes through a 
‘maximum. For large angles, y >>h/jl, the extensional energy becomes so 
: large that »x¢ >>, so that » approaches ¥,. 

The transitional interval is centered around an angle which is roughly 

proportional to h/l. As h tends to zero the shell becomes inextensible, but it 

retains its bending modes, if in the limit v9 is kept constant. 

Figure 8 clearly illustrates the effect of the competition between extensional 
and inextensional displacements for minimizing the total potential energy. 
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Zusammenfassung 


Im ersten Teil (I) werden Formeln entwickelt, welche die kinetische und poten- 
tielle Energie diinner Schalen beliebiger Form explizit als Funktion eines belie- 
bigen Verschiebungsfeldes ausdriicken. Dies geschieht auf rein analytische Weise 
ohne Zuhilfenahme geometrischer Darstellungen einzelner Schalenelemente. Von 
besonderem Wert fiir viele Anwendungen ist es, dass vollstandige Freiheit in der 
Wahl der Flachenkoordinaten besteht. 

Im zweiten Teil (II) wird zundchst die Anwendung der Theorie auf zwei 
bekannte Beispiele gezeigt: 1. Platte in schiefwinkligen Koordinaten, 2. zylin- 
_drische Schale. Als neues Beispiel wird 3. die schraubenflachige Schale eingehend 
 behandelt. Fiir einen besonderen Einspannungsfall wird die Frequenz der Grund- 
schwingung als Funktion der Verwindung angendhert numerisch berechnet. 
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On the Theory of Disc-Loaded Waveguide 


By HEINRICH DERFLER, Ziirich!) 


§ 1. Introduction 


An electromagnetic difference-equation of importance in the theory of 
traveling-wave tubes was established in a preceding paper [1]?). The parameters 


of this equation were derived from a variational principle resulting from the - 


Maxwell-Lorentz field theory. Periodical structures of cylindrical symmetry 


were considered in general. It is the purpose of this paper to apply those 


methods to a special structure and to calculate explicitely the related para- 


meters. The theory of the disc-loaded waveguide was considered before in the 


literature [2], [3], but little care was taken to the effect of wall thickness. Since — 


this is closely related to the problem of maximum energy interchange with the 
electron beam it will be considered in detail here. In addition the concept of 
coupling impedance will be criticized. The notations of [1] will be used here and 
(1)-(24) refer to equations used in that paper. 


§ 2. The Evaluation of the Admittance Function 


The admittance function G(6) was defined by 


+D/2 +D/2 
BZ, 2440) = (Z EZ adh az 
2, _ IO Se 9 
G(6) £5 a pA _ ~Pl2 -D/2 


G +D/2 +D/2 : (15) 
i) i} E*(Z ) dZ dz’ 
—D/2 —D/2 


The kernel related to the disc-loaded waveguide is by (6) and (9): 


TAS = On ome | (Mn A) Y, AC B) Pe rAGs B) Nas A) 
BAe 2) eee Dy tea Oe Bree a ee ee 


X COS (F Als 5) COSH 7% (= + 5] (25) 
+00 


1 Ji(% A) 17) 2nna+ 0 
or ap ee ee 


n= — 


1) Institut fiir Hochfrequenztechnik der ETH. Paper supported by a research fellowship of the 
Battelle Memorial Institute. 


*) Numbers in brackets refer to References, page 206. 
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Since the problem was shown to be stationary with respect to small variations 
in E(Z) a suitable approximation will be taken as E(Z). Such an estimation is 
suggested by electrostatics in the following form: 
2Z\2)-Cs D 

E De G+ GC, l1- (+) | IZ|<=. (26) 
C, should be chosen such that E(Z) is correct near the edges R= A, Z = +D/2 
of the slots. It may be shown by the methods of conformal mapping [4, page 325] 
that near a right angular corner C, = 1/3. This value would lead us to Bessel 
functions of the fractional order 1/6. Since tables of these functions are not 
available at present we choose C; = 1/2 which is the correct value near a 


knife edge. Hence ealata ie 
E(Z) = E(0) [1- (45) |" 121 Aah (27) 


We have to evaluate the following integrals [5, page 151] 


+D/2 


[az E@) COSN (5 ~ 5) =F (0) en Joma) (-1)", n=2m, 
- D2 
= Q; n=2m+1, (28) 
+ D/2 


‘ 2 +06 Du 2nn-+ 60 
[az BW) expiZ **F—— = E00) oe I(t bl. 


~ Dye 
Introducing these and (25) into (15) gives: 
BA \ 20  Ji(B A) ¥o(B B) — Jo(B B) Vi(B A) 
C BD Jo(B A) ¥.(B B) — Jo(B B) Yo(B A) 
Je(n mt) Ji(Hen A) HY(Man B) — Jo(Hon B) H{)(Han A) 
Poe a en DT A) HiP\aee BY = Jltan BY HOM an A 


G(6) = 


2nnu+ 0 
ene Ym L Jo(% A) Ja 2L D) 


(29) 
- Now we make the practical assumptions: 


po= 24 <an, © <az. (30) 


Since } > a this yields the following approximations 
aaa pA /1- (aa) epee Rt Si p>, 


=pBy/i-( Gey etn 0% 


Han B 


Han D=BD ed n> 0. 


¥ 
| 
| 
y 
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From the asymptotic representations of Bessel functions [5, page 137] we have ; 


n>0O 


Ji(ton A) HY (Hen B) — Jo(Hon B) HO (Mon A) 
Jo(Han A) AY (Men B) — Jo(Hon B) HY (Hen A) b 4 < 22 
—a 


I 
: 


hence: 


4n . Je(n 4) : Si(Han A) AY (Usn B) i Jo(Hon B) HWY (Un A) 
n=1 Hon D Jo( Han A) AY (Mon B) y= Jo( Hen B) AW (Hen A) 


(31) 
5? Je(m 2) . 
Se | 
‘Using the abbreviations 
2 6 22 ad 
e=BpA=—", z= >5—,- a= 7, b= 48% (32) 
and introducing 
v4 EY oof 
mL = L|/p2— ( ‘ns ) ‘a 2n/(Z) = (n+ 2)? 
as well as the approximation (31) into (29) gives 
_% | ma f(x) You dja) — Jo(x b/a) ¥,(x) y7 Ja (n 2) 
WES Gat Thaler bla) = hla Balhae, = gee 
+00 Jule aig ait at | 
peme Veley® — (wn ya Jolm + Pee 
Finally the admittance function may be written in the form 
G(0) == = R(x) — S(x, 2) (33) 
where by definition 
= Se 
_— 1 Sil*) Yo(* b/a) — fy(x b/a) Y,(x) 
RO) = ila) Vola Ba) = Jule Bla) Yala)” 34) 
yr Sala )/Jolo 1a ae 
) § (2 2) 
(x, 2) = 2. ae = Join + 2 Ck Goa la (35) 


n=1 


a ee 
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The function R(x) is readily computed from tables [5], but the series in S(x, z) 


converge rather slow. A series transformation will be applied to represent the 
result in a form which is more convenient for computation. 


§ 3. The Transformation of S(x, z) 


KUMMER’S transformation of infinite series will now be applied to S(x, z). 


f ‘ . . . . . 
_ The idea is to substract a series the terms of which are asymptotically propor- 


tional to those of S(x, z) and to add again their sum. Using the asymptotic 
representation of the Bessel Functions [5, page 137] we have 


NAC: RAC V ) il ee 0(n-2) 


V(%/a)?— (n+ 2)? In +2! 


Je(n + z 8) = a |n| >co. 


|n| >oo, 


_ Hence (35) will be written in the form 


si. — 32 | OVUM) 1 _ | ag +38) | 


p Hn | Vx)? — (m2? Ie +2 | Jo 


In fact the terms of the first series decrease as 0(n~%) which is suitable for com- 
putation. If the sum of the last two series were known, the problem would be 
solved. Again these series converge rather slow and another transformation is 
necessary to prepare them for computation. This will now be done by means 
of the Mellin-inversion formulae. From [7, page 198] we have 
k+ioo 
f(a. = oS : 2-1 7(5) Ries) c=" (== =| wens. Os wee by. 
k-100 


Replacing s > s — 1 and x > x & gives 


k+100 : 
r(x OD 1 ~ 2h\an0 1 s—1 cep SN 
ee oz | GG) ite Te-9r (AGA) was 
sit Le weee,, 
taking x = n + z and summing over yields 
pres Jin +2 9) 
n=0 n+ Z 
k+t0o Ez = 
ee fi (5) 5 r(- : )re s) I-(=5 =) tis, 2)ds 1<k<2. 
k—100 (37) 
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To save writing we put the integrand of (37) 


1s) = 5 (5) 1S re-9 r-(- 


This is a meromorph function the poles of which are located at 
ie td LAr A irae AES 


The corresponding residues are 


Res I(s) = In (5) ey ety | 


=U 
11 (9\2" I(2n) P(2n + 1) NGS | 
Res Ls) = (—1) : i I (n+ 1) an(Z | Dai ue 
—1\7 Bale Aa 
Res I(s) = ( 2 ‘ey 1 ee C(2.n, 2h a> 0: 


These will be used to evaluate (37) by means of CAucHyY’s theorem. 
Extending the path of integration to the left half plane yields (Figure 1) 


k+100 32/2 


1 ies 3 
rai | Ts) ds=Res 1s) + 3’ Res I(s Oey Coa e'* dy, 


2N—1<r<2N41. 


It can be shown, using STIRLING’s formula for Js) [8, page 279] and HERMITE’s 
formula for ¢(s,.z) [8, page 270] that the last term vanishes if N co, Hence 
by (39), (38) and (37) 


\n- an I'(2n) I(2n+ 1) (40) 
+31 ; (5 ) rae ye? P,n(z) 


Oesy =e 


Replacing z > 1 — z and adding to (40) gives 


fi Zz P) 
Dy fee) : = 21n (5) — 27 - pe—1) — y(-2) 


= n-1(8\2n F(2n) (2 1 41 
+23, (-ye+(F)" Teg Pent) pig f 4M 


*) t(s, z). (38) 
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Now from StirLine’s formula and from | P,,,(z)| < 4/(2 2)?” 


= 2, (=)*" 
a ne \ 5% : 


| op T(2n) F(2n+1 


) 
2 I4(n + 1) Prn(2) 


Thus (41) converges if # <z and the new series is usefull for computation if 


O< 2. 


ims 


Figure 1 
The singularities of I(s). 


I = path of integration in (37); J7 = contour yielding a convergent series; 
III = contour yielding an asymptotic series. 


We may extend the path of integration to the right half-plane also (Figure 1) 


k+ic + 2/2 

1 i 

aa | 1 — SP Res Is) + ae [terete dp, 2N<r<2N42. 
k—-ioo uae 


Again using the formulae of STIRLING and Hermite it can be shown that for 
any fixed value N both Res I(s) and the last integral can be made arbitrarily 


small by making # sufficiently large. Hence by (39), (38), and (37) the following 
asymptotic series representation is obtained: 
+ 


om Jé(m+z8) 1 3 (—1)"-1 (2) C(2m,2). (42) 


N+2 28 ht (2 — 1) 
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Espécially for 2:10 = a 


yr Je(n an) 1 1 ! Pa Mites 2 4 
Poe: ~ 6 720 * 4480 sel Se 
Finally, using (36), (41), and (43) we have 
2 1 
S(x, 2) = 2In() — 2y + > — ve -1) - y(-2) 
= nai fo\2n (2 n) (22 + 1) 
+ 2 OUR Cera i hee 
yt | Sle VV Jole V_) Lal |, pares 
' zo 
ee Galt mae tay {0 
Onl (0S = x 
which is fairly good for computation. 
§ 4. The Optimum Values of the Slot Parameter 
In our notations the Pierce gain-parameter reads [9, page 50} 
te, ait 
C.= 7 (K Ayre: (45) 


Hence C3 is a quantity which should be as large as possible. From model 
considerations similar to those used in [9, page 36] the following proportionality 
may be shown to hold: 


D 
Z = const 3 (46) 


the coupling impedance Z being essentially the same for all the space-harmonic 
waves 0) + 2a (24). 

For an electron beam filling the drift-space we have the transadmittance 
[9, page 23] 


Kearedila) tte ay | By eens gael Ps i se 


which is essentially the same for both klystrons and traveling-wave tubes. As 
a consequence of the assumption E(Z) = E(0), |Z|< D/2 the factor 
[(sinD/2)/(D/2)}2 appears in the quoted paper in place of J?(D/2). To be 
consistent with the principle of stationarity (13) these calculations were 


“Vi a 


id 
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repeated yet using (27) in order to prove (47). Hence for a given space- harmonic 
wave in synchronism with the electron beam + L=0,+2nz [9, page 37] 
the factor (D/2) J>(D/2) should be as large as possible, yielding the condition 
(D/2) J,(D/2) = (1/2) Jo(D/2), D/2 = 0-9408; 3-9594; ... The corresponding opti- 


- mum values of the slot parameter are 


Gan4 — 77 0.0408, 47 3.9504,... 9<a. (48) 
Usually the optimum should coincide with the middle of the band-pass response 
6, = 2/2 giving L = | 4 + 1 | a/2 and the values of the slot parameter # listed 
in Table 1. 
Table 1 
Optimum values of the slot parameter }. 


n | 0 =nall 
—2 0-538 2-263 
—1 1-254 -- 
1 0-753 _ 
+ 2 0-418 1-760 


Small values of # strongly reduce the coupling between adjacent slot-fields 
yielding small bandwidthes or very high electron velocities to satisfy the broad- 
band condition [9, page 43]: electron velocity = phase velocity = group velocity 


L=0+2nn=05% n>0. (49) 
Hence, in computing S(x, z) from (44) 0 = 1 and # = 2 were used giving nearly 
optimal values of the gain parameter C, for m= —- 1 respectively n = + 2. 
By this procedure the structures obtained from the tables in § 7 should be 
suitable equally well to both, forward- and backward-wave interaction. 


§ 5. The Evaluation of the Electronic Interchange Parameters 


The phase-shift per elementary cell can be infered from (24) G(9) = 0 which 
by (33) gives 
Ze R(x) = S(#, 2). (50) 


This equation may be solved graphically as shown in Figure oe 
An example useful for broad-band amplification by interaction with the second 
forward-space-harmonic wave n = + 2 is shown in Figure 3; the corresponding 


quantities are listed in Table 2. 
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Figure 2 


Evaluation of bandpass-response « = 10, # = 2. 


Figure 3 


The characteristic quantities of the disc-loaded waveguide. 


a = 10; = 2; b/a = 2; synchronous velocity: B = 1/15-17;n = + 2. 
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Table 2 
The characteristic quantities of the disc-loaded waveguide. «a = 10; } = 2; bla = 2. 


Having found the phase shift per elementary cell, the coupling impedance is 
readily obtained from (24) Z = 2/G'(6)) which by (33) is 


ve 4a€ ec) (51) 


x Oz 


This quantity can be computed by use of the tables § 7. The results are shown 
in Figure 3, the values being listed in Table 2. 
The disturbing impedance is given by (24). 


cotg es + dW 


wa i 0, sin Oy 
ON Pa G’(0,) [cots 2 cos — cos | : 


- From (33), (50), and (51) we deduce 


W=rara 2% smart 12 cotgnz + OW. (52) 
The complex roots z, = 9;/22, k > 0 of (50) and the corresponding derivatives 
0S/dz should be known in order to evaluate dW. The mathematical situation 
may be seen from Figure 4. 

It seems to be hopeless to find the z,° at any reasonable expense. They 
correspond to sets of space-harmonic waves which are below cut-off in the range 
of frequencies under consideration x < /}. 
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/Slx,z)/ 


ip 


Figure 4 
Relief of S(z, z) = | S(x, z) | e¥ (*. 2) 
for a fixed value of x < p, (diagrammatically) « 29° pes VP —p?, Tog eee 


A very rough approximation of their contribution to W may be infered from 
asymptotic formulae using «z,~(k+ 1/4)ia, k>O and neglecting slow 
varying quantities. The result will be given here without proof: 


2 23 2 2 ; 
jie One o (5) e~ (52/24) (7+ 8) | (53) 


% Sale Am Pe oa 


The contribution of (53) to (52) is shown in Figure 3, the values being listed in 
Table 2. Now the Pierce space-charge parameter can be computed from [9, page 


28] yielding ore 
Once — [sin? (32) : 1 3 “| ; (54) 


Fortunately the first term coming from space-charge usually predominates and 
hence the error from (53) may not be serious. 


§ 6. The Concept of Coupling Impedance 


Form formal considerations, using the theory of three-terminal-pair junc- 
tions [9, page 26] the coupling impedance is suspected to be 


Z,= 1b. (55) 


7 
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; where by definition app 
V = [EQ dz = E00) &, (56) 
—d/2 
SE \« 
PaahiZie arf aif ze Ha2nurdr, (57) 
ede 0 


Using the decomposition into space-harmonic waves, the field components 


_ were shown to be [1]. 


+D/2 se 
pert Hy ’ ' 1 Si(%m = re, {ie 2nn+@ 
Ho = E(Z\ dZ e2. She OMY ORD UC 2) 
~ D2 
1 0H » 
a] Ne 


F . oe 2 0 . flv, R) Serr 2nu-+ 0 
IEA oe as tk: aoe exp i(Z’ — Z) 7 . 


—D/2 AES 


Now using (27) and (28): 


He= 7h. (0) 
<3 te fee n(22* 0) ex (12 285*?), 
E, = 4" £(0) 
ae iim ) p(2% 5+" p) exp | an 


From Poyntinc’s Theorem we have the power-flow in z-direction 


es +00 
1 
Pi) = 7 [Et Ho2nrdr=— 3) Pam(Z) 
0 i 
where 
dn\2 Il(2na+ CS ENE BE dal 
Pn, m(Z) = af ca. ay) ETO. Ay alms £ 


(58) 


A 
x EXP Bat (ma MZ f Jalvn R) lm R) RAR + 0 
0 
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but 
nL 2 Berke 2 +@ 
1 nN I 
Poo. + aad Sa eae 
-t2 
(59) 
‘ Silt A) — Jol A) Jo(%m A) (et 6 D) = Mm: 
Tilt A) ak 
==), n=m. 


This indicates a strong interaction of the space-harmonic modes within an 


elementary cell when a perturbation is present. The interection is such that the — 


continuity of fields across the boundary R = A is assured, resulting into a 
definite ratio of power carried in the modes: 


Pun _ th 2nn+0 Ji(2na+ 6 D/2L) 
mm me 2mat+O 722mn+ 6 D/2L) 
Tilt A) = Jol%m A) Jo(%m A) J8(%m A) 
Ji(Ym A) = Jol%m A) Jo(%m A)  J8(% A) * 


Hence, if one of these modes interacts by synchronism with the electron beam, 
all the others are produced as well; the total power-flow being: 


(60) 
x 


P= PZ) = 3) P= a Arh (2)? ExQQ | 
yr 2nat 6 Jim A) — Solr A) Alo 2 6 a 
WG i(% na+ 
x 2, iD Tiga een | 
Introducing this and (56) into (55) gives 
1 _ B 
Zz, mz A® ca | 
oo + (G2 
x 2m et 0 S84 A) ~ Hmm A) ln A) pf 2nate py | © 
poco ne Ji(% A) Jal 20 ) 
Now the correct value, predicted by consequent field theory is from (29): 
es OGe) 
fo LO 
B Ty 2nnx+0 J2v_ A) — Jolrm A) Im A) B/2na+ 0 
= plaid Me as RUN ZY Jo\ln 2\'n 
aA er BLS Tits, Jo aa D) 


gf ®) oat 2nu+ 6 2n2z+ 06 
ae ae a Jol om fs D) hi Pius D) 


. 


= 
f 
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and by comparison with (62) 


1 1 


p $s 


where 


Pe eD pS A Til A) 2nn+ 0 g +6 
a aed Mm L Iol% A) | iE D) A <i aes ). ee 


Hence, a shunt impedance Z, appears across the formally defined coupling 
impedance Z,. This was to be expected because the theory of three-terminal- 
_ pair junctions assumes interactions between nearest neighbours only [9, page 27]. 


—co 


§ 7. Tables of the Functions R(x) and S(x, z) 


The functions R(x) (34) and S(x, z) (44) were computed from tables [5] 
whereas 0S(x, z)/0z was infered from S(x,z) by numerical derivation. The 
error is less than + 5 units in the last place. 


10 R(x) 
En bla 
2nalA 1-8 1-9 | 2-0 | Q1 | 2.2 | 23 Q-4 
1-00 5-805 
105 4-637 
1-10 4-684 3-602 
1-15 4-833 3-747 2-670 
1-20 3-978 2-904 1-824 
1-25 4-286 3-208 2-138 1:042 
1°30 4-672 3°574 2-508 1-435 0-3140 
1-35 5-144 4-009 2-930 1-868 0-7790 
1:40 4-513 3-405 2-340 1:274 0-1660 
1-45 5-107 3-940 2-854 1-798 0-7232 
1-50 4-554 3-418 2-348 1-294 0-2023 
1:55 4-049 2-940 1-883 0-8233 
1-60 3-588 2-501 1-451 0:3806 
1-65 3-164 2-093 1-047 
1-70 2-773 1-716 0-6658 
1-75 2-411 1-362 0-3075 
1-80 2-075 1-031 
1°85 1-760 0-7172 
1:90 1-467 
1-95 1-189 
2-00 0:9281 
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Seay) ea Aare — 0), i sare — Il 


x=2nalA 


0-00 6°743 7-498 7-938 8-582 9-578 11-261 14-643 
0-05 6:593 7:284 7-679 8-252 9-123 10-547 13-242 
0-10 6-210 6-749 7-019 7-465 8-074 8-977 10-565 
0-15 5-730 6-111 6°312 6-585 6-962 7:508 8-304 
0-20 5-265 5-524 5-655 5-827 6-054 6-358 6-777 
0-25 4°872 5-047 Briley 3242 5-381 5-561 5-795 
0-30 4-563 4-689 4-743 4-816 4-907 5-020 5-161 
0-35 4-330 TAZo 4-466 4-518 » 4-581 4-657 4-750 
0-40 4-180 Zo 4-282 4-322 4-371 4-428 4-495 
0-45 4-091 4:151 4-178 4-211 4-252 4-301 4-355 
0-50 4-061 4-119 4-144 4-175 4-215 4-259 4-309 


= 0S(#; 2)/02 «=2na/l= 10, P= ndii=1 


zg x= 2nalA 
6/22 0-0 | 1-0 1-2 | 1-4 | 1-6 | 1-8 | 2-0 
0:00 0-000 .- 0-000 00-00 00-00 00-00 00-00 00-00 
0-05 5-729 8-152 10:22 12-37 16-90 26-48 48-65 
0-10 9-076 12-40 14-34 17-69 22-99 32-04 52-15 
0-15 9-723 12:58 13-98 16-74 20-55 26:52 37-48 
0:20 8:672 10:72 11-88 13-44 15-74 19-20 24-33 
0-25 7-026 8-296 9-062 10-01 11-30 13-08 15-64 
0:30 5-331 6-179 6-588 7-141 7-865 8-847 10-15 
0:35 3-805 4-338 4-556 4-869 5:264 5-794 6-488 
0:40 2:427 2:705 2°841 3-025 3°237 3-482 3-853 
0:45 1-179 1-303 1-360 1-448 1-531 1-662 1-824 
0-50 0-000 0-000 0-000 0-000 0-000 0-000 0-000 


S(#,2) «=2aa/l1=10, ®=and/l=2 


x= 2nala 
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0-00 0-000 
0-05 6:355 
0-10 10-12 

0-15 10-89 

0-20 9-776 
0-25 8-010 
0-30 6-130 
0-35 4-403 
0-40 2-820 
0-45 1-376 
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—OS(%, 2/02 a= 2 24/2 = 10, o= nd/2 = 2 


x=2malaA 


| 1-2 | 1-4 | 


00-000 00-00 00-00 


00-00 00-00 
11-05 13*23 17-91 27-70 50-28 
5252, 19-00 24°41 33313 53-94 
15-30 18-08 21-91 27-88 38°89 
13-05 14-62 16-89 20-29 Phaye Al 
10-02 10:95 12:29 3-99 16-47 
7358 7-902 8-607 9-552 10-81 
3129 D445 5°821 6-326 6-964 
3:22) 3-398 55996 3-848 4-178 
1-558 1-621 1-720 1-835 1-990 
0-000 0-000 0-000 
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List of Symbols 


Only symbols additional to those used in [1] are listed. 


By, 

Cy 

K 

P,(2) = (2+ B)"/P(m + 1) 
Q» 
*a=BPpA=2a2a/A 
Z)= Iy/Vo 

we 0/2 7 

“a= 2ra/l 

Ts) 

y = 0:577216 

C(s, 2) i; am (n ae Fat 
o= 2d/ ‘ 

A 

yp(z)= = In (e+ 


Bernoullian numbers [6, page 245] 
Pierce gain parameter 
Transadmittance [9, page 23] 
Bernoullian polynomials [6, page 541] 
Pierce space-charge parameter 
Frequency parameter 

Beam impedance 

Phase-shift parameter 

Drift-space parameter 

Evuter’s Gamma-function [8, page 235] 
EuLeErR’s constant [8, page 235] 


Generalized Zeta-function of RrEMANN [8, page 265] 


Slot parameter 
Wave length in free space 


Gauss’ Psi-function [5, page 16] 
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Zusammenfassung 


In einer vorangehenden Arbeit [1] wurde eine elektromagnetische Differen- 
zengleichung aufgestellt. Die auftretenden Parameter folgen aus einem Varia- 
tionsprinzip, das aus den Maxwell-Lorentzschen Gleichungen gewonnen wurde. 
In vorliegender Arbeit ist diese Methode auf den periodisch mit Lochblenden 
belasteten Wellenleiter angewandt. Die Ergebnisse beanspruchen praktisches 
Interesse im unteren Millimeter-Wellenbereich. 
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Thermostatische Zustandsainderungen des trockenen und des 
| nassen Dampfes 


Von Lane S. DzunG, Baden?) 


1. Einleitung 


Die Eigenschaften des nassen Wasserdampfes fiir massige Driicke im ther- 

modynamischen Gleichgewichtszustand konnen heute wohl als gut bekannt 

_ betrachtet werden, nachdem eine Anzahl Dampftafeln und Enthalpie-Entropie- 
Diagramme mit iibereinstimmenden Zahlwerten in diesem Bereich schon vor- 
handen sind. Die Zustandswerte sind gewoéhnlich in einer verwickelten Weise 
miteinander gebunden, so dass es bis jetzt nicht méglich war, die Gesetzmassig- 
keiten der Zustandsanderungen durch einfache Gleichungen zu formulieren. 
Diese Schwierigkeit spielt zwar fiir die einzelnen Berechnungen keine Rolle, 
denn man wird lediglich ein Zustandsdiagramm zur Hand nehmen und alle 
nétigen Zahlen ablesen. Fiir die allgemeinen Untersuchungen wird dies jedoch 
oft unangenehm. Die Durchrechnung aller Méglichkeiten wird sehr zeitraubend. 

- Ein Problem kann nicht mehr in seiner Allgemeinheit unabhangig vom Zu- 
standsort gelést werden. 

Eine der Schwierigkeiten liegt darin, dass in die thermische und die kalo- 
rische Zustandsgleichung sowohl die Eigenschaften der fliissigen und der gas- 
férmigen Phase als auch die transzendente Dampfdruckformel eingehen. Wir 
werden in dieser Arbeit zeigen, dass man nach einem neuen Verfahren zur 
Berechnung der Zustandsinderungen diese Schwierigkeit leicht umgehen kann. 
Durch geeignete Definition einiger neuer Zustandsvariablen und Differential- 

_ koeffizienten ist es méglich, auch fiir den Nassdampf einfache Gleichungen zu 
entwickeln, die denjenigen fiir ideale Gase ganz ahnlich sind. 

Weitaus schwerwiegender ist das Auftreten der Unterkithlung des Dampfes 
und der Uberhitzung der Fliissigkeit, das heisst Ubersattigung. Die Ursachen 
solcher metastabiler Zustande sind zweifacher Art, naimlich die Oberflachen- 
spannung und die Relaxationszeit. 

Die Beitrage der Oberflachenphase zu den thermodynamischen Grdssen 
miissen beriicksichtigt werden, wenn die Fliissigkeit in sehr viele feine Tropfen 
verteilt ist. Der Dampfdruck, oder, genauer, die Fugazitat bei kleinen Tropfen 
ist grésser als bei grossen Tropfen um einen Faktor, welcher exponentiell mit 
der Oberflachenkrimmung des Tropfens zunimmt. Dieser Ubersittigungsfaktor 
wurde schon von Stopota [1]?) eingehend untersucht. Die Rechnung zeigte, 


1) Aktiengesellschaft Brown, Boveri & Cie. 
2) Die Ziffern in eckigen Klammern beziehen sich auf das Literaturverzeichnis auf Seite 222. 


dass der Einfluss dann bemerkbar wird, wenn der Tropfenhalbmesser in die 
Nahe des Hundertfachen eines «Molekiilhalbmessers» (~ 2 A) und darunter 
kommt. Da die kleinen Tropfen sich mit der Zeit zu grésseren Tropfen zusam- 
menschliessen, ist dieser Zustand metastabil. 

Diese mikroskopische Erklarung der Metastabilitat kann man durch eine 
phanomenologische erganzen. Die Einstellung des thermodynamischen Gleich- 
gewichtes durch Phasenumwandlung, darunter auch das Verschwinden der 
Oberflachenphase bei der Bildung grosser Tropfen, bedarf einer gewissen 
«Relaxationszeit ». Fassen wir den Umwandlungsprozess als eine Art chemische 
Reaktion auf, so kénnen wir eine Verdampfungslaufzahl x definieren. Sie ist 
nichts anders als die bekannte Dampftrockenheit. ° 

Bei zu rascher Zustandsanderung, zum Beispiel beim Durchgang durch eine 
Stosswelle, bleibt x praktisch konstant, ahnlich wie die Reaktionslaufzahl 
einer gehemmten Reaktion. Die Partialpotentiale der gasformigen und der 
fliissigen Phase sind verschieden. Die Differenz ist eine Art Affinitat, wovon 
die Umwandlungsgeschwindigkeit :abhangen diirfte. 

Das hier berichtete Rechnungsverfahren des Nassdampfes ist allerdings nur 
fiir den thermodynamischen Gleichgewichtszustand giiltig mit der Gleichheits- 
bedingung des Partialpotentials der fliissigen und der gasférmigen Phase. Tat- 
sachlich gibt es in der technischen Anwendung Faille, bei denen diese Gleichheit 
besteht, zum Beispiel wenn der ZeitmaBstab der Zustandsanderung geniigend 
gross ist gegentiber der Relaxationszeit. Bei der Bestimmung der Maschinen- 
grosse nach wirtschaftlicher Betrachtung bedeutet die Variation der Zustands- 
grossen keine wirkliche Zustandsanderung, sondern nur einen rein thermostati- 
schen Vergleich benachbarter Zustande. Hier ist die Gleichgewichtsbedingung 
sicher erfiillt. Bei der Expansion des Dampfes in den Niederdruckstufen einer 
Turbine ist der ZeitmaBstab der Zustandsainderung jedenfalls nicht sehr klein 
gegentiber der Relaxationszeit. Die Schaufelerosion bestatigt das Auftreten 
von Kondensation, das heisst Einstellung des Gleichgewichts innerhalb dieser 
Zeitdauer. Jedenfalls kann man bei der Berechnung der Zustandsdnderung in 
Turbinen den Gleichgewichtszustand als einen idealen Grenzfall betrachten. 
Denn jede Verzégerung in der Einstellung des thermodynamischen Gleichge- 
wichtes bringt immer einen Entropiezuwachs oder einen Verlust an nutzbarer 
Energie durch den nachtriaglichen Ausgleich der Partialpotentiale. In diesem 
Sinne ist auch das mit einem Zustandsdiagramm graphisch bestimmte « Warme- 
gefalle» im Nassdampfgebiet zu verstehen. Der wirkliche Verlust wegen Poten- 
tialausgleiches ist allerdings durch die «Wasserbremsung» verschmiert. Wir 
werden am Schluss diese Ausgleichsverluste noch etwas naher betrachten und 
ihre Grdésse nach der Methode der Thermodynamik irreversibler Prozesse ab- 
schatzen. 

Im Laufe der Entwicklung zeigt sich, dass wir zu diesem Zwecke einige all- 
gemeine Beziehungen der Thermodynamik mit Vorteil etwas anders als bisher 


i 
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formulieren kénnen, worauf die Anwendung auf den Sonderfall des nassen 
Dampfes miihelos und unmittelbar folgt. 


2. Der isentrope und der isenthalpe Exponent 


Wir betrachten ganz allgemein ein geschlossenes System von einer Stoff- 
_komponente mit zwei Freiheitsgraden. Ein solches Gebilde nennen wir einen 
allgemeinen Dampf. Es kann sich dabei sowohl um die rein gasférmige Phase 
handeln, in diesem Fall sprechen wir von einem trockenen Dampf oder einem 
realen Gas, oder um einen nassen Dampf, das heisst die Gleichgewichtsmischung 
- der fliissigen und der gasférmigen Phase. Im letzten Fall kommt das Massen- 
verhiltnis der beiden Phasen, ausgedriickt durch die Dampftrockenheit, als ein 
neuer Freiheitsgrad hinzu; dagegen wird wegen der Gleichgewichtsbedingung 
die Temperatur durch den Druck gegeben. Man hat daher in jedem Fall zwei 
Freiheitsgrade. 
Wir verwenden den Druck # und das spezifische Volumen v als die zwei 
Zustandsvariablen und definieren in bekannter Weise den isentropen Expo- 


nenten dlogp 
= 0§ 
en nee it ; () 


Dann gilt fiir kleine Zustandsanderungen langs einer Isentrope offenbar 
p v* = const (s = const) . (2) 


Ganz ahnlich definieren wir den isenthalpen Exponenten 
bas dlog p 
erste gee ), ; 8) 
womit fiir kleine Zustandsainderungen langs einer Isenthalpe gilt: 


po™=const (=const). (4) 


Es lasst sich nun zeigen, dass bei irgendeiner Zustandsanderung das Ver- 
halten des Systems durch k und m eindeutig beschrieben wird. Diese Exponen- 
ten sind Zustandsfunktionen und kénnen im allgemeinen Fall von zwei Zu- 
standsvariablen abhangen. 

Wir definieren eine Zustandsgrésse 

cate (5) 
durch die Gaskonstante R. Fiir ideale Gase ist t offenbar identisch mit der 
Temperatur 7. Fiir Nassdampf ist die Temperatur T schon durch den Druck p 
bestimmt, 7 ist aber eine unabhiangige Zustandsgrosse. 
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Aus diesen Definitionsgleichungen kann man nun ganz allgemein die Be- 


ziehungen * 
vdp=dh=—~—— Rdt (ds 2=0) (6) 


fiir eine Isentrope und | 
vdp=—-Tds=—™"_ Rd (dh =0) (7)% 


fiir eine Isenthalpe ableiten. Man muss hierzu beispielsweise nur die Gleichung _ 
(1) in der Form 


vdp+kpdv=0 (ds=0) 
schreiben und Gleichung (5) sowie die Beziehung 


dh = Tds+vdp (8) 


beriicksichtigen. Gleichung (6) ist iibrigens schon in der elementaren Thermo- 
dynamik bekannt. Wir schreiben die Gleichungen (6) und (7) noch in einer 
etwas praziseren Form: 


oh Os 
(ie), =48) T(5e),--#R ©) 
mit den Abkiirzungen 
Pieri: WE 
a ei fen rey (10) 
Nun kénnen wir t durch / und s ausdriicken: 
Ot OT 
ar (on), 4 + (5), @ 
und mit Riicksicht auf (9) die Beziehung 
1 1 
Rat me dh — - Las (11) 
erhalten. 
Fiir die technische Anwendung schreiben wir zuerst rein formell 
dh=gpvdp, —Tds=(1—g)vdp. (12) 


Die hier definierte Grésse y kennzeichnet die Richtung einer Zustandsinde- 
rung. Die beiden Gleichungen sind wegen (8) aquivalent. Die Grésse v dp in 
Gleichung (12) ist offenbar die innere quasistatische Durchstromarbeit. Das 
Differential dh ist die 4ussere Durchstromarbeit im adiabaten Fall. » ist daher 
der Wirkungsgrad bei einer infinitesimalen Expansion oder die Reziproke des- 
selben bei Kompression. Wir definieren nun die polytrope Zustandsanderung 
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~ durch die Bedingung @ = const und den polytropen Exponenten durch 


age (sce? Le (13) 


- Dann gilt fiir kleine Anderungen 
pv" =const (gy =const) . (14) 
- Mit Gleichung (5) erhalt man 


vdp=———, Rat (y = const) (15) 


_ wie in den Gleichungen (6) und (7). Durch Einsetzen von (12) in (11) hat man 


aber 
1 


je ae) v dp 


R dt = ( 


und durch Vergleich mit Gleichung (15) und wegen (10) mit » = /(m — 1): 


gag 1—@ 1 7 1—@ 
- Fs WAIT oder a ate (16) 


m 


und in etwas gelaufigerer Form: 


n—1 Ne il m — 1 
Hy b + oS a. 


nN 


Die Zusatzbedingung gy = const ist hier nicht mehr notig, da  selber in der 
Gleichung vorkommt und somit die Richtung der infinitesimalen Anderung 
vorschreibt. 

Fiir den spiteren Gebrauch leiten wir noch eine Beziehung ab durch Elimi- 
nation von dh zwischen Gleichung (8) und (11) und Division durch t nach 
Gleichung (5): 
ae 2) ds he RP —RG-1) ae Re te 


C 


Wir werden diese Gleichung bald benutzen kénnen. 
Bis jetzt sind unsere Ausfiihrungen vollig allgemein. Es sind keinerlei Gas- 
eigenschaften angenommen worden. 


3. Das ideale Gas 


Fiir das ideale Gas gilt definitionsgemass p v = R T oder t = T. Die thermo- 
dynamische Identitat 
(Sp), = — Ger) 
Op! Tr OT] p 
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wird in diesem Fall ie , 
eae “gig 
Nach Gleichung (17) haben wir aber die allgemeine Beziehung 
(55) we) R/p ee oe ee 
ph ~ ~ (TH 1 Ap) i ae S| 


Damit die beiden Ausdriicke fiir t= T identisch sind, muss aber = co und 
wegen Gleichung (10) m = 1 sein, wenn die andere Lésung A = 0 wegen innerer 
Instabilitat fiir endliche Temperatur ausgeschlossen wird. Aus Gleichung (9) 
hat man dann, wegen , 


oh) _ _ (dh/dr), _ AT _ 
cas (ds/dt), , 
dh=RdT, A=AKT)=—?, h=A(T). (18) 


Die zweite und die dritte Gleichung folgen daraus, dass in der ersten die linke 
Seite ein vollstandiges Differential steht, weshalb die rechte Seite integrierbar 
ist. c, ist die spezifische Warme nach der iiblichen Definition. 
Aus den Gleichungen (16) und (17) hat man ferner 
n—1 k—1 1 


- Me a 
n ok eae oa oodles 


aT ad 
a ee (19) 


Die Entropie ist separabel. Diese Beziehungen sind aus der elementaren 
Thermodynamik bekannt. 


4. Der ideale Dampf 
Wir definieren den idealen Dampf durch die Beziehung 
m=1, w=oco. (20) 
Dann gilt nach Gleichung (9), wie bei (18), 
ah=ARdt, A=Xt), h=h(z). (21) 


Ferner gilt nach Gleichung (16) und (17): 


Mat Fail: 1 2) 
aot apt, 1a8, 2) 


T dt dp 
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Wegen der zweiten von Gleichung (21) ist die rechte Seite von (23) separabel 
und integrierbar. Die linke Seite ist daher ein vollstandiges Differential, was 
bedeutet, dass der « Kompressibilitatsfaktor» 


_ nur von s abhangig sein kann, also 
2 = 2(s). (24) 


Definiert man nun wegen (24) eine Grésse o durch do = ds/z, dann wird aus 
Gleichung (23) 


dt d 
do = RA“* —R = (25) 

Die Gleichungen (21), (22) und (25) sind denjenigen fiir ideale Gase, Glei- 
chungen (18) und (19), véllig 4hnlich, wenn T durch 7, s durch o und c, durch 
AR ersetzt werden. 

Es ist méglich, umgekehrt aus Gleichung (24) allein alle andern Gleichungen 
dieses Abschnittes zu gewinnen. Die Beweisfiihrung ist etwas verwickelter. 
Wir schreiben Gleichung (11) in der Form: 


Care ee raga 


und erhalten durch kreuzweise Differentiation 


R (Se) = Lae (Ge). (26) 
Einsetzen von Gleichung (24) ergibt: 
R (2) = = [PO + Sa (27) 


Anderseits haben wir aus Gleichung (17): 


1 


et ee 
T z 


ne 
Die rechte Seite ist bei Giiltigkeit von Gleichung (24) sofort integrierbar. Damit 
ist die linke Seite auch ein vollstandiges Differential. Die Bedingung hierfiir 


ist aber R [2 (=)], = [se (sa). 


oder nochmals mit (24) 


bah Jee eR 
7 z 


R eas pe Eek 


~ \ 0s = dt 
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Ein Vergleich mit Gleichung (27) liefert sofort die Bedingung 


oder wieder Gleichung (20), wenn wir den Fall A = 0 wegen innerer Instabilitat 
ausschliessen. 

Dieser Beweis fiir das Aquivalent der Gleichungen (20) und (24) ist allge- 
mein giiltig. Fir trockenen Dampf kann etwas anders vorgegangen werden [3]. 
Fiir den Sonderfall, wenn ausser Gleichung (24) noch k= const gilt, hat 
TRAUPEL [4] die Ahnlichkeit eines solchen polytropen Dampfes mit einem 
polytropen Gas, dasheisst einem idealen Gas mit konstanter spezifischer Warme, 
zuerst erkannt und auf den Umstand hingewiesen, dass verschiedene technisch 
wichtige Dampfe (Wasserdampf, Luft) tatsachlich diese Eigenschaften in guter 
Annaherung besitzen?). 

Bei konstantem k und 4 kann man (18) und (21) integrieren und erhalt 


h=ART + const (28) 
fiir das polytrope Gas und 
h=4 Rt+const (29) 


fiir den polytropen Dampf. 

Die Analogie zwischen einem idealen Gas und einem idealen Dampf darf 
nicht zu weit getrieben werden. TRAUPEL [4] hat schon darauf hingewiesen, 
dass das Aquivalent dann aufhért, wenn Warmeleitung?) ins Spiel kommt. 
Der Grund ist, dass durch Gleichung (25), in der Form dh = t do + v dp ge- 
schrieben, zwar sowohl dem ersten Hauptsatz als auch einem Teil des zweiten 
Hauptsatzes der Thermodynamik Geniige getan wird, jedoch nicht dem nullten 
Hauptsatz [2]. Denn zwei Systeme im Gleichgewicht haben zwar dieselbe 
Temperatur 7, aber nicht immer dasselbe r3). Aus demselben Grund kommt 
auch der Grosse o nicht dieselbe Bedeutung zu wie der Entropie s. 


5. Der semiideale Dampf 


Der Schritt vom idealen Gas zum idealen Dampf kann nochmals fortgesetzt 
werden, indem wir statt (20) oder (24) die weniger strengen Bedingungen 


A=i(t), w= p(s) (30) 


1) Unser idealer Dampf ist etwas anders als derjenige nach Lrtp [6]. Der polytrope Dampf 
ist jedoch identisch mit seinem idealen Dampf bei konstanter Oszillationsenergie. 

2) Warmekonvektion stort das Aquivalent nicht, wohl aber Warmestrahlung. 

8) Im Sinne des Theorems von CARATHKODORY bedeutet dies folgendes: t ist zwar ein integrie- 
render Nenner fiir ein homogenes System mit zwei Freiheitsgraden, jedoch nicht fiir ein System 
mit mehr als zwei Freiheitsgraden, zum Beispiel fiir mehrere Systeme im diabaten Gleichgewicht. 
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auferlegen. Die rechte Seite von (17) wird nun separabel und die linke Seite 
vollstandig. Damit ist der Koeffizient von ds nur von s abhangig oder 


tend orx cecil ke beara ga 


Wegen des Verschwindens der linken Seite von (26) ist das Produkt A T/u nur 
von s abhangig. Die letzte Gleichung wird 


1 ($2) eee ib Oz\ Az 
22 NOT \c. fk tT? aoe 


_ Wegen der zweiten von (30) besteht diese Abhangigkeit auch fiir das Produkt 
AT =At/z, solange (20) nicht gilt: 


lee al dlc eae 
Elimination von (0z/0r), liefert mit (30) 


di dt Aaa. 
ee 
i, t lb T 


Die linke Seite ist wieder separabel und die rechte Seite vollstandig. Damit 
kann 4/u oder nochmals mit der ersten von (30) auch w nur von t abhangen, 
jw = w(t). Das ist aber nur bei konstantem jw mit der zweiten von (30) vertrag- 
lich. Zusammen mit der letzten Differentialgleichung erhalt man dann die Be- 
dingungen aa 
ana G coust, 7 = const. (31) 
Durch (31) wird der Funktionenbereich von (30) eingeengt. 

Wir nennen ein gasférmiges Medium einen «semiidealen» Dampf, wenn es 
nur Gleichung (30), jedoch nicht Gleichung (20) oder Gleichung (24) gehorcht. 
Die praktische Anwendung ist besonders dort zu finden, wenn die Exponenten 
Aund w annahernd?) konstant sind, also wenn es sich um einen «semipolytropen » 
Dampf handelt. Dann lassen sich viele Ausdriicke explizit integrieren. Die ver- 
~ mehrte Komplikation gegeniiber dem polytropen Dampf ist unbedeutend. Wir 
werden sehen, dass dies fiir den nassen Wasserdampf zutrifft. 

Wir erhalten aus den Gleichungen (2) und (6) fiir die Isentrope identisch 
wie bei einem polytropen Dampf: 


p, vt = py 05, Pies (=y', ies (2) (s = const) ; (32) 


Pe To V4 T, 


ia Rie ia AR, (Se) —1] (s =const) (33) 


1) Der Fall, dass A und y genau konstant sind, ist nur bei «4 = co mit der Integrierbarkeit von 
Gleichung (17) vertraglich. 
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und aus den Gleichungen (4) und (7) fiir die Isenthalpe: 


pw =p,0%, t= ( “4 Ves w2 (= le (a= const). (34) 


P2 Te: Vy T2 


Das Entropieglied lasst sich nicht explizit integrieren. 
Fir die Polytrope gilt nach den Gleichungen (12), (14) und (15) 


ptt=p%, B=(B), 22=(2)™ (@=const); (35) 


P2 T2 Vy T2 


ha—-hky=QvR(t,.—-t) =pvRy lea - 1] (py =const) (36) 


1/>. 


mit 


Wir eliminieren ds zwischen den Gleichungen (8) und (11) und erhalten statt — 


Gleichung (17) 


1 1 v 
Rév=(z—=) dh+— ap. (37) 
Daraus ergibt sich fiir die Isobare 
He il Ory T, 
R(ty—1)=(7- |) a—h), G2==* (=const). (38) 


Der polytrope Dampf ist in diesen Gleichungen durch die Spezialisierung 
auf m = 1, 4 = oo gegeben. Dabei werden die Gleichungen (33), (36) und (38) 
in bezug auf die Beziehung zwischen / und t identisch. Man kann dann nach 
Gleichung (29) die Enthalpie # durch 7 allein ausdriicken, was bei einem semi- 
polytropen Dampf nicht méglich ist. 

Die Grundgleichungen der Str6mungsmechanik haben dieselbe Form fiir 
polytropes Gas und polytropen Dampf [4], bis auf den nullten Hauptsatz. Fiir 
einen semipolytropen Dampf versagt dieses allgemeinere Aquivalent, da Glei- 
chung (29) nicht mehr gilt. Nur wenn die Zustandsainderungen bestimmten 
Gesetzen folgen, sind die entsprechenden Formeln, Gleichungen (32) bis (38), 
anwendbar. In der Technik kommen aber die meisten Probleme unter diesen 
Kategorien vor. Fiir solche Probleme sind die angegebenen Gleichungen von 
grossem Nutzen. 


6. Der trockene Dampf 


Fiir den trockenen Dampf oder das reale Gas gelten die Gleichungen in 
Abschnitt 2. Definieren wir: 


a & = (See 
ee 2 abe inge ib B 


Ill 


ale Pe ie (582), = a (39) 
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dann kénnen die folgenden Beziehungen abgeleitet werden: 


‘ne ee Be Ue clas We pa Le 1 
VN ree ee, EF cael =). (40) 


-Daraus ergibt sich aus Gleichung (16). Zum Beweis leitet man am einfachsten 

zuerst den Ausdruck fiir 1 ab [5]) und erhalt k und m durch Spezialisierung 
-auf py = 1 und p = 0. 

Aus (37) und (5) und mit den thermodynamischen Identitaéten findet man 


ifaae) = R (55), (Gr), (30), (or), = oS) ee 
: und mit (40) die friiher [5] abgeleitete Beziehung 


r (op). ~ p(t 3): 


7. Der nasse Dampf 


Fiir den Nassdampf verlieren alle Grdssen der Gleichung (39) den Sinn, 
da p und T voneinander abhangig sind. Um die Ausdriicke fiir k und m zu 
finden, beschranken wir uns auf solche massige Driicke, die geniigend weit vom 
kritischen Druck entfernt sind. Insbesondere betrachten wir immer die gas- 
férmige Phase als ein ideales Gas und die fliissige Phase als isovolumetrisch. 
Die Zustandswerte der fliissigen und der gasformigen Phase werden allgemein 
durch die (oberen) Indizes f und g bezeichnet. Der Strich und der Doppelstrich 
bezeichnen iiblicherweise die entsprechenden Phasen im Sattigungszustand, 
das heisst x = 0 und x = 1. 

Mit der Abkiirzung Av = v" — v’ schreiben wir das spezifische Volumen des 
Nassdampfes in der bekannten Form*) 


gene’ +(l— xv =v +2 do 


und wegen der Annahme v’ = const 


dv Ox dv” 
—-) = — —_. 42 
(sp), 4° (ap). + * ap 42) 
Fiir den Gleichgewichtszustand mit zwei Freiheitsgraden ist « durch # und s 
gegeben. Man kann daher schreiben 


( Ox (0s/0P) 
Op 

1) Dort [5] wird yp = 1/pstatt p verwendet. 

2) Die additive Eigenschaft aller extensiven Grissen fiir Nassdampf ist exakt, solange die 
Oberflachenphase nicht beriicksichtigt wird. Dann gibt es auch keine Wechselwirkungen zwischen 
den beiden Phasen. 


aes (Os/Ox 


e 
p 
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Wegen 


hat man aber 


( 


as 
Ox 


), $< 6 = As nnd (s5),> x 


Nach unserer Annahme iiber die Eigenschaften der beiden Phasen haben wir 
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s=xs"+(1—2)s' 


as’ 
dp © 


ds” = 
ap + (L—=4) 


unter Beriicksichtigung der Clausius-Clapeyronschen Gleichung: 


ds”, 0s8 Os9 ate ies 4 lp Av 
‘dp (or lnee inch dp i i As 2 
ds’ _ dst! aT oo’ dw 
dp dt ip Fr is 
dv” __(0v9 ($F) ah q! he Me roe : Av 
ae OT)» dp p Aes 


Der Index g bei der spezifischen Warme c, ist weggelassen, da keine Verwechs- 
lungsméglichkeit besteht. Einsetzen in (42) ergibt 


OO \SE | ety Spat Cy Av Ae i of. An eee 
eat see ae Pas) ~ O— (e) + 2 ( As me 
und nach Gleichung (1) mit pv” = RT: 
a aa R Beye Ag hie a5 (Ao)? 0. chy Tene 
Rk v [ae a ant (=) yt R |+ vu” R (G-) < (43) 
Den Ausdruck fiir m erhalt man am einfachsten iiber Gleichung (41): 
i ils il PORTA Cp Ain RO Bee a eae (44) 
wien es ned po Peele eas oF Asset Ase ) 
Fir v’ <v" kann man Av durch v” ersetzen. Somit wird 
Ling ean” R R\2 c’— cy a ele Ne bed R 


Pir v' << xv" kann man den Faktor x v"/v streichen und v” /v durch 1/x ersetzen. 
An der unteren Grenzkurve!) bleibt nur das letzte Glied von (43): 


Jenseits der Grenzkurve wird nach unseren Voraussetzungen k = 
x=O—, und k= k’, m=m! nach Gleichung (40) 


k 


R 


A ela DA ASHEAS 
R c 


u 2 


(x = 0+). 


durch die Grenzkurven machen k und m begreiflicherweise einen Sprung. 


1) Hier miisste natiirlich die Beschrankun 


den, was grundsatzlich keine Schwierigkeit bereitet. 


(46) 


m = oo fiir 
fiir x= 1+. Beim Ubergang 


g auf isovolumetrische Flissigkeit aufgehoben wer- 
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Q 


Fiir Wasserdampf im Atmospharendruck und darunter sind die Werte von 
k und m nach Gleichung (45) berechnet und in Kurvenform (Figuren 1 und 2) 
angegeben. Es ist auffallend, dass die Werte fiir « > 0,6 in diesem Druck- 
bereich erstaunlich konstant sind. Sie kénnen angenahert als 


Le wat fie fe 
ce 2 


Zeuner ee abode 
as ee ee oem na 
jes mae 


SS 
| VAN 
V/ 


0 0 7 
0,5 x 1,0 0,5 x 1,0 
Figur 1 Figur 2 
Isentroper Exponent des nassen Isenthalper Exponent des nassen 
Wasserdampfes. Wasserdampfes. 


angenommen werden. Der nasse Wasserdampf verhilt sich in diesem Bereich 
wie ein semipolytroper Dampf. Dieser Zustandsbereich ist aber derjenige, wie 
er in den Niederdruckstufen der Kondensationsturbinen und in Kondensatoren 
vorkommt. Fiir 0,4 < « < 0,6 wird sogar m = 1, das heisst, er verhalt sich wie 
ein polytroper Dampf mit k » 14/13. Gleichung (24) ist fiir trockenen Wasser- 
dampf sehr gut erfiillt, jedoch nicht fiir nassen Wasserdampf?). 

In Figur 1 sind auch die Werte nach der von STODOLA zitierten, von ZEUNER 


stammenden Gleichung 
k-=1,035 + 0,1 x 


aufgetragen. Diese empirische Gleichung ist fiir den Bereich 0,4 er ee 
geniigend genau. Sie versagt jedoch vollstandig fiir kleinere Werte von *. 

Die Gréssen 2 und p weisen selbstverstandlich eine grossere Veranderlich- 
keit auf als & und m. Indessen kann gezeigt werden, dass‘es fiir die Berech- 


1) Vel. hierzu Figur 2 von TRAUPEL [4]. 


220 Lane S. DzuneG ;  ZAMP” 


nungsergebnisse meistens nur auf k und m oder auch w/A ankommt. Die erste 
Bedingung von (30) ist fiir das Gebiet x > 0,6 sehr gut erfiillt. Auch die erste 
von (31) ist mit hoher Genauigkeit befriedigt, wenn man fiir ~ nicht den jewei- 


ligen Wert, sondern eine geeignet gewahlte Konstante, etwa den Wert fps bei | 
x = 0,8, einsetzt, was wegen der zweiten Bedingung von (31) vollkommen be- ; 
rechtigt ist. Der Ausdruck 1 A/(uo, — A) betrégt namlich iiberall 150 K mit — 


Schwankung von + 10 K innerhalb des genannten Gebietes. Diese Priifung 


ist eine weitere Stiitze fiir die « Semiidealitat » des Nassdampfes. 
Die Schallgeschwindigkeit a im trockenen Dampf ist bekanntlich durch 


C=RPT=SER TT (47) 
gegeben. v und R sind hier auf die Masseneinheit bezogen. Im nassen Dampf 


gilt dies nur dann, wenn die Schallperiode gegeniiber der «Relaxationszeit» 
der Phasenumwandlung sehr gross ist, was fiir die hérbaren Schallwellen nicht 


zutrifft. Trotzdem kann man auch hier mit Vorteil eine Geschwindigkeit nach — 


Gleichung (47) definieren, die wir die Pseudo-Schallgeschwindigkeit nennen. 
In gewissen Anwendungen kommt dieser Geschwindigkeit auch dieselbe Be- 
deutung zu wie der Schallgeschwindigkeit im Trockendampf. Sie ist eine be- 
queme Bezugsgrésse, und man kann mit ihr auch eine Pseudo-Mach-Zahl bilden 


wie mit der wahren Schallgeschwindigkeit. Diese letztere wird vielmehr gegeben — 


durch a* = (0p/0@),,; ohne Gleichgewichtsbedingung, oder 
ao Gln. = or Gel] =e Goh 


Gash p vy" fir of << xp 


und 


8. Ubersattigung 


Bis jetzt haben wir nur den thermodynamischen Gleichgewichtsfall unter- 
sucht. Dem Nassdampfsystem mit zwei Phasen wurde die Zusatzbedingung 
uw’ = w (Gleichheit der Partialpotentiale der fliissigen und der gasférmigen 
Phase) auferlegt und damit ein Freiheitsgrad entzogen. Wir werden nun diese 
Zusatzbedingung aufheben und nehmen nur thermomechanisches Gleichgewicht 
an. Der Druck und die Temperatur sind zwar fiir die beiden Phasen gleich, 
jedoch nicht durch die Dampfdruckformel gebunden, 

Fiir unser abgeschlossenes, das heisst geschlossenes und adiabates, im 
thermomechanischen Gleichgewicht stehendes System kann sich die Entropie 
nur durch Umwandlung zwischen den Phasen andern. Die Gibbssche Funda- 
mentalgleichung wird einfach 


lds. A dx. (48) 


: 


: 


x 
" 
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A ist hier die Affinitat der Verdampfung und ist gegeben durch die Differenz 
des Partialpotentials der beiden Phasen, A = — (wu? — w’). 
Wahrend der Verschiebung des Verdampfungsgrades oder der Dampftrok- 
_kenheit x von einem beliebigen Anfangswert auf den Gleichgewichtswert x° 
bleibt das System energetisch abgeschlossen, das heisst, es wird weder Warme 
noch Arbeit mit Aussen ausgetauscht. Die Energie und das Volumen oder, fiir 
‘den stationaren Fliessprozess, die En- 
thalpie und der Druck+) bleiben un- 
_verandert. Die Temperatur des Systems 
andert sich vom Anfangswert JT auf 
den Gleichgewichtswert T°. Das ist 
die Sattigungstemperatur entsprechend 
dem konstanten Systemdruck #, bezo- 
gen auf sehr grosse Tropfen. 
Es ist méglich, die Abhangigkeit der 
_ Affinitat A vom Verdampfungsgrad x 
unter diesen Bedingungen eindeutig 
auszudriicken. Die durch T dividierte 
rechte Seite von (48) lasst sich dann 
mit gewissen Annahmen iiber die spezi- 
fischen Warmen elementar integrieren 
und ergibt die Entropieanderung beim 
Einstellen des Gleichgewichtes. Der end- 
_ giiltige Ausdruck ist etwas kompliziert. Figura 
In gewissen Fallen gelangt man nach Adiabate Expansion des Nassdampfes mit 
der Methode der letzten Abschnitte und ohne Ubersattigung. 
schneller zum Ziel. Wir suchen die En- 
tropiezunahme bei der isentropen Expansion des trocken gesattigten Dampfes, 
Zustand 1, zuerst ohne Phasenumwandlung, auf einen gegebenen Enddruck, 
Zustand 2 und mit nachtriglicher (isenthalper) Kondensation auf den Gleich- 
gewichtszustand 3, bezogen auf sehr grosse Tropfen. Figur 3 zeigt diese Zustands- 
- Anderung auf einem modifizierten Mollier-Diagramm [7], in welchem eine Isobare 
2S des Trockendampfes ausser derjenigen des Nassdampfes 3 S eingezeichnet ist. 
Punkt 4 ist der Endzustand bei quasistatischer Expansion auf denselben End- 
druck, das heisst der Endzustand nach dem iiblichen Dampfdiagramm. 
Nach Gleichung (33) kann man sofort schreiben: 


ly — hy =A Rt, [1 - ey") | I — hy =I Ro, [1 - (oy) 


1) Der Einwand von BINNIE und Woops [8] gegen die isobare Annahme hat hier keine Giultig- 
keit; da wir nur die Stagnationszustande betrachten und damit das Problem auf dasjenige ver- 
nachlassigbarer kinetischer Energie reduzieren. 


222 . Lane S. DzunG . ; ZAMP 


A und 4° sind die entsprechenden Exponenten fiir trockenen und fiir nassen 
Dampf. Da die Nassdampfisobare gradlinig ist, hat man mit der Sattigungs- 
temperatur 7° beim Enddruck 3 = fy: 


T° (sg — Sg) = T° (Sg — $4) = hy — hg = hg — Ig 


d 

Penman eb ("af") | 
Ao Lenght Tay APL EDS i 
ee k ko py ae 


Das ist der Verlust an Verfiigbarkeit bei Einstellung des Gleichgewichtes und — 


gibt gleichzeitig den Entropiezuwachs an. 


9. Zusammenfassung 


Es wurde die Berechnungsmethode der thermostatischen Zustandsande- 
rungen in einem allgemeinen dampfformigen System nochmals tiberpriift. Durch 
Einfthrung des isenthalpen Exponenten neben dem tiblichen isentropen Expo- 
nenten kann die Rechnungsmethode in einer einheitlichen Form, ahnlich wie die 
bekannten Formeln fiir ideale Gase, gebracht werden. Der Sonderfall des poly- 
tropen Dampfes und des semipolytropen Dampfes werden besprochen, wobei die 
aus der elementaren Thermodynamik bekannten Formeln fiir polytrope Gase 
mit unwesentlichen Anderungen sich anwenden lassen. Die Ausdriicke dieser 
Exponenten sind fiir den trockenen Dampf und den nassen Dampf angegeben. 

Damit haben wir ein analytisches Rechnungsverfahren sowohl fiir den 
trockenen als auch fiir den nassen Dampf. Auf eine Niitzlichkeit des analy- 
tischen Verfahrens wurde schon frither [4] [5] hingewiesen. 

Das Problem der Einstellung des Gleichgewichtes wird kurz gestreift. 
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Summary 


The change of states of a closed thermodynamic system having two degrees 
of freedom can be expressed analytically by means of the isentropic and the 


_ isenthalpic exponents, suitably defined. By imposing certain plausible conditions 
on these exponents, it can be shown that many well-known formulae, originally 


derived for an ideal gas with constant specific heat, become applicable more 


generally. Low pressure wet steam of not too low dryness fraction as occurred in 


the low-pressure stages of condensing turbines obeys these imposed conditions 
fairly well. 
The problem of supersaturation and of the drift toward equilibrium has been 


_ briefly touched upon. 


(Eingegangen: 26. Mai 1954.) 


Two Dimensional Flow with Constant Shear Past Cylinders 
with Various Cross Sections 


By Anprew R. MitcHett and JAmeEs D. Murray, St. Andrews, Scotland?) 


Introduction 


The stream function y for the rotational flow of an incompressible inviscid 
fluid in two dimensions satisfies the equation 


ory | Op et 
Ox2 Oy? +o=0, (1) 
where the vorticity is given by 
Fy 94x 
‘Ox Oy” (2) 


gz and q, being the velocity components parallel to the x and y axes respectively. 
Methods of solution of (1) have been given by TsIeEN [1]*), Coomss [2], and 
James [3] for flow with constant vorticity past stationary cylinders with circu- 


" lar and airfoil cross sections. The forces on the cylinders have been calculated 


using the methods described but except in the case of the circle, closed expres- 
sions for the stream function in the physical plane can be obtained only with 
great difficulty, if at all. Accordingly, in the present paper it is intended to 
obtain the stream functions for constant shear flow past a variety of cylinders 
using in each case an extension of TsIEN’S method based on the natural co- 
ordinates of the system. From each stream function, the main characteristics 
of the constant shear flow can be determined. 


1) The University, Department of Mathematics. 
2) Numbers in brackets refer to the References, page 234. 
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As far as the present authors are aware no theoretical solutions have yet 
been obtained for flow with variable rotation past cylinders in two dimensions, — 
although several investigations using approximate methods have been carried 
out. These include an evaluation using finite difference methods by one of the _ 
present authors [4] of the rotational field downstream of the bow shock wave | 
formed when a parallel supersonic flow impinges on a blunt nosed two dimen- © 
sional obstacle. | 


The General Transformation and the Circular Cylinder 


An auxiliary stream function Y which satisfies LAPLACE’s Equation is 
introduced by means of the transformation 


ax?t by? (3) 


il 
F=ypt zo atb Z 


where a and 0 are constants. For the case of flow with constant shear parallel 
to the x-axis 


il 
ap WE Merrainign ONES (4) 
where U is the velocity along the x-axis, and from the above transformation 


a (x? — y?) 


1 
T= U. paiee@ ES a (5) 


The stream function for shear flow past a circular cylinder derived by © 
TstEN can be obtained very quickly by putting a= in (3) and (5). This 
simplification gives Y for parallel flow as the imaginary part of the complex 
function U z+ 1 w 22/4, where z= x +14, which according to MILNE THom-, 
son’s Circle Theorem [5] is modified to U (z + 1/z) + (i w/4) (22 — 1/z?), when 
a circle of unit radius is inserted at the origin. Substituting Y the imaginary 
part of this latter complex function into the transformation formula, the requi- 
red stream function 


(6) 


i = = w+ U sind (r - =| _ _ 7) (2 v2 sin? + er) 


is obtained for rotational flow past a circular cylinder, the stream function 
taking the value zero on the cylinder. 

For all cross sections examined in this paper, other than the circle, it is 
convenient to put a = 0 and so from (3), 


1 
ytleedt aia Ee (7) 


where Y’ is a solution of LAPLAcr’s Equation. 


> - 
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The Elliptic Cylinder and the Flat Plate 
(A) Derivation of the Stream Functions 


| In the case of flow with constant shear past an ellipse centred at the origin 
_ (Figure 1), elliptic co-ordinates £, are introduced by means of the formulae 


x=ccosh€cosy, y=csinhésinn , (8) 


Figure 1 


Illustration of rotational flow past an elliptic cylinder. 


where c is a constant. The ellipse & = & (major and minor axes a and 6 along 
the x- and y-axes respectively) is taken as the solid boundary and the velocity 
profile at x = — oo is given by (4). Expressing LaApPLAce’s Equation in elliptic 
co-ordinates and solving by separation of variables the function VY satisfying 
the boundary conditions at infinity is found to be 


w—Uy+ JS (A, cosny + B,sinnn)e""® (n=0,1,2,...) 9) 
n=0 


which on substitution in (7) gives the stream function p. In order that y may 


retain the constant value zero on the ellipse § = & O0SnS2 gt), the coeffi- 


ZAMP VI/15 
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cients in (9) must take the values 


A,=A,=Ay=4g="-=0, Ag=4+ 508%, Ag= ae ea 
Bio ue(stn\", EE HE oe iy. 
leading to the required expression 
y= - wo b? + U[csinhé — be™-*] siny 
(10) 


nie - w[c? sinh2é 4-{ 63 e2Go-s) _ c®'sinh?é } cos2 n| 


for the stream function. This value of y satisfies the three necessary conditions © 


(1) pro =0, 
(2) w=0 on €=6,, 
fl 1 
(3) Yr Oe Ui ar OF as E>oo. 


In an exactly similar manner, if the velocity profile at y = + co is obtained 
from 


p=Vi- Fox, (11) 
the stream function for the shear flow past the ellipse € = &, is found to be 


y= Toa + Vie coshé — a e»-@) cosy | 
(12) 
Serr [c? cosh2é — { az e7(So-) _ 62 cosh®é }cos2 n|, | 


where again y takes the value zero on the cylinder. In the limit as the minor 
axis tends to zero, the ellipse becomes a flat plate of length 2 c, and the stream 
function for the shear flow past this plate is given by 


y= 2 wc?+ Vecsinhé cosy | 
(13) 
iGo c*[cosh2é++ (sinh2 & — sinh?) cos2 7]. 


The flow fields represented by (10) and (12) are symmetrical about the y 


and x axes respectively. The velocity components in the € and » directions are 
qe and g,, where 
Q=—o SS, = iat Se Aa 14) 
a g(§ 7) On In g(é,7) o0& (14) 
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pad g(é, ) = ¢ (sinh2é + sin?7)". 


The stagnation points are given by 


(B) The Elliptic Cylinder with its Major Axis in the Direction of the Stream 


In the case of the ellipse, whose major axis is in the direction of the stream, 
there are two stagnation points in the range 7 <7 S27 for all values of the 
non dimensional parameter (U/w 6). Denoting this number by N, the positions 
_ of the stagnation points, which always lie on the ellipse, are given by 


(15) 


sing = 7 N-— = (N?+ 2b ne 


GD 


- where 0 < N < oo. The range of movement of the points is from positions on 
the ellipse, the joins of which to the origin make angles 
; b 1/2 
fuel pas =) 
with the major axis, to the ends of the major axis. The depth of either stagna- 
tion point below the x-axis, using (15), is given by 


= 5 NZ (Nt+y)" (16) 


where 6’ = y/b and y = 2 b/(a + 6) lies in the range 0 <y <1. The graph of 
6’ against 1/N for y = 1/2 is shown in Figure 2. The corresponding result for a 
circular cylinder is obtained by putting y = 1 in (16). The non dimensional 
ordinate of the stagnation stream line at « = — ce is given by 
y= N— (NP + 5) 
where y’ = y/b. 
The number of stagnation points in the range 0 = 7 < a depends on N. If 


2a 4-6 
ON Seto’ 


there are two stagnation points in the range. They lie on the ellipse and their 
positions are given by 

: i) 1 2 26 \1/2 

Sit) = cx N+ z (v a =m 5] 


It 
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Figure 2 


Variation of the position of the lower stagnation point for the circle, ellipse, and flat plate. 


there is one stagnation point which lies on the y-axis and whose position is 
given by 6 < y < oo. The stagnation points thus start off from positions on the 
ellipse, the joins of which to the origin make angles 


' b 1/2 
oo Gove 

with the major axis, move round the ellipse away from the major axis until they 
coincide at the upper end of the minor axis, whence the combined point moves 
off along the y-axis to infinity. The flow field when the upper stagnation points 
have just coincided is shown in Figure 3. In view of the symmetry of the field 
about the y-axis, it is sufficient to illustrate the region x < 0. 

The resultant force on the cylinder acts in the y-direction, since the flow 
field is symmietrical about the y-axis. Its magnitude per unit length of cylinder 


is given by + 0/2 


Lae ao/ a _¢, sinn dn, 


— 7/2 


where g:_, the value of the velocity on the ellipse, is obtained from (14), and 
g is the density of the fluid. This integral can be evaluated to give 


L==¢,0 (a+b) U?, (17) 


where the lift coefficient c, = 2 a/N, with N = U/(w 6) as before. 
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Figure 3 
Illustration of the stream lines when the upper stagnation points coincide. 


(C) The Elliptic Cylinder with tts Minor Axis in the Direction of the Stream 


An exactly similar state of affairs exists for the ellipse whose minor axis 1S 
in the direction of the stream. Starting with (12), formulae are derived as before 
with a and } interchanged and V replacing U. The flow is now of course sym- 
metrical about the x-axis. The value of N at which two of the stagnation points 
coalesce is (a + 2 b)/2 (a + b) which reduces to the limiting values of 3/4 and 
1/2 for the circular cylinder and flat plate respectively. The depth 6’ of the 
lower stagnation points below the y-axis (x <0) is again given by (16) with 
N=Vjao, 6’ =x/a, and y = 1, 2a/(a + 6) (1 <y< 2), and 2 for the circle, 
ellipse, and flat plate respectively. The graphs of 6’ against 1/N for the above 
cylinders (y = 3/2 for the ellipse) are shown in Figure 2. It should be noted that 

_ 6’ and N change their significance at y = 1 in this figure. The lift force on the 
cylinder, which acts in the x-direction is given by (17) with N = V/a o. 


The Parabolic Cylinder 


(A) Derivation of the Stream Function 


In the case of constant shear flow past a parabolic cylinder (Figure 4), 
parabolic co-ordinates w, v are introduced by means of the formulae 


x=ur—-v?, y=2uv. (18) 
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Figure 4 


“Illustration of rotational flow past a parabolic cylinder. 


The parabola v = v, is taken as the solid boundary and the velocity profile at 
x = — oo is again given by (4). Proceeding exactly as in the case of the elliptic 
cylinder, the required stream function is found to be 


yp=2UUuv—20u? v? + [{a(n) cosun + b(n) sinunbe-*"dn, (19) 
0 


where m the separation constant is continuous and takes all values between 
zero and infinity. In order that y may take the value zero on the parabola 
v= Up, the coefficients A(n) and B(n) must satisfy the equation, 


co 


[{A@ cosun + B(n) sinu nhdn =2wyui—-2U yu, 
0 


where A(n) = a(n) e~”” and B(n) = b(n) e~*”". Using Fourier Transforms, 
A(n) and B(n) are determined as 


oO [o.0) 


An) = = @ “8 | u*cosnudu, B(n) = — 2. U v | w sinnudu. (20) 


0 0 


In order to keep the coefficients finite, the integration in (20) is restricted to a 


sl elie 


«4 
~ 
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; finite range A which can be chosen arbitrarily large. With this limitation, (20) 
_can be integrated by parts to give 


j 4 1 : 

} A(n) = — ov [= (a2— =) sind n+ <> cosa], 

. (21) 
B(n) = — Uv ee sind n — Ze cosd n| , 


| Substituting these values in (19) and after considerable manipulation using 
DE HAANn’s Tables of Integrals [6], the required stream function is obtained as 


| y= 2Uur—20u2o2— 9 M—) (7 _ 2 v4) 


n 
(vu — U9)? + (u — A)? 2u 
ee er ete ale 
a 2A (uv — Up) 4 w v2 A (v — Up) 
> SE eae ee aL an a pee 


This value of y satisfies the three conditions 


(1) V2yt+o=0. 
(2) yw=0 on v= providedu <A. 
(3) pruy-Zoy as Uv —>Ooo. 


In verifying conditions (2) and (3), it should be kept in mind that 


2A (v — %) 
(v — v9)? + u® — A? 


te 


takes the values z and 0 for v = v9 (wu <A) and v = on respectively. 


(B) Examination of the Flow Field 


The stagnation points, which lie on the parabolic cylinder v = v9 (u<A) 
below the x-axis are given by 


y , r / 1 re U 
(N + 2’ v9) [20 vg + (1— uw") wplog——i +m (l—w *)|—202=0, (23) 


where N = U/(w A2), u’ = u/A (<1), and vj = v/A. There is one value of w’ 
satisfying (23) for each N in the range 0 S N Soo. The stagnation point at 
N = 0 occupies a position on the parabola depending on vj, and as N increases 
it moves along the parabola towards the x-axis reaching the latter at infinite NV. 
The non dimensional depth of this stagnation point 6’, is given by the value wv’ 


i 
| 


~ - j 
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| 
which is a solution of (23), since 6’ = 6/(2 A v9) = u’. The graph of 6’ against Nis _ 
shown in Figure 5. | 
The stagnation points which lie on the cylinder above the x-axis are given by — | 
ee ae 
t= 


(N — 2’ vf) [2 w’ of + (1 — w'?) 09 log tau’ (L—w'?)|+2052=0. (24) } 


-O5 


-04 


Figure 5 


Variation of the position of the lower stagnation point for a parabolic cylinder. 


If 0 < N< vw, there is one stagnation point on the parabola above the x-axis. 
When N = 0, this stagnation point and the one below the axis are symmetri- 
cally placed with respect to y = 0. As N increases, the upper point moves along 
the parabola away from the x-axis. Its position at N = vj depends on vj. When 
vy = 1, the stagnation point moves from a position w’ = 0-39 to wu’ = 0-80 as 
N goes from 0 to 1. For N = vj, the situation is rather more complicated. For 
example when vj = 1, at N = 1 a second stagnation point appears above the 
x-axis at uw’ = 1. As N goes from 1 to 1-15, the original point moves from the 
position wv’ = 0-80 to u’ = 0-95, whilst the new stagnation point moves from 
wu’ =1towu' = 0-95. The points thus coalesce at uw’ = 0-95 and for 1:15< N < oo, 
there are no stagnation points on the cylinder above the x-axis. 


(C) Discussion of Results 


The results obtained in this section require careful consideration in view of 
the presence of singular points on the parabolic cylinder at u =A. A guide to 
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_ the accuracy of the method can be obtained by deriving the flow field past a flat 
_ plate by the methods of this section and comparing it with the field obtained 
_ by considering the plate to be the limiting case of an elliptic cylinder. 
Solving LAPLACE’s equation in cartesian co-ordinates by separating the 
variables, and proceeding as before, the stream function 
1 x x2 + (yi — A)? 
p=Uy- Foy U-oy) log Po | 
(25) 
1 2A % woAx 
ie. aye = Se 8 <i 
[Uy 7 oly x?)| tg arte yee a 


is obtained for shear flow past a flat plate which lies along the y-axis in the 
range —A<y<+A, where / again is an arbitrarily large finite quantity. The 
points y = + Aare singular points of the field. 
The main features of the two flat plate flows are indentical. The lower 
stagnation point in each case starts off at a position on the plate at N = 0 and 
as N increases the point moves along the plate towards the axis reaching the 
centre of the plate at N = oo. The stagnation points in each case are symmetri- 
cally situated on the plate at N = 0 and as N increases the upper stagnation 
point moves away from the axis leaving the plate at N = 1/2. For 1/2 <NS es, 
the top half of the plate is free from stagnation points. The values of the 
velocity at identical points in the flow fields are not in so close agreement. 
Sufficient to say that the closest correspondence occurs on the plates away 
from the edges. 
As a result of this comparison, it seems reasonable to suppose that the flow 
field represented by the stream function (22) will contain all the principal 
features of the field of flow past a finite parabolic cylinder (v = v),0 Su SA). 
In addition, (22) can be used to give a good approximation to the flow on the 
surface of the parabola in the vicinity of the nose. As A can be chosen arbitrarily 
large, the valid region can be extended considerably. 


Conclusions 


It is clear from the foregoing study that the power of TstEN’s method for 
obtaining constant shear flow patterns past cylinders is greatly increased by 
using if possible the natural co-ordinates of the system. 

The use of elliptic co-ordinates for example enables stream functions to be 
obtained for flow past elliptic cylinders (excluding the circular cylinder) and the 
flat plate. The elliptic co-ordinates, however, are of no advantage in solving 
the problem of shear flow past one branch of a hyperbolic cylinder. The slit 
along the x-axis between the foci of the families of hyperbolae divides each 
branch of every hyperbola into two sections each with its own value of 7. For 
example the branch of the hyperbola passing through the point c/2 on the 
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x-axis has equations y = 2/3 and y = 5 2/3 above and below the x-axis respec-_ 
tively. As a result, it is not possible to obtain the required stream function | 
using the methods outlined in the present paper. | 

Using parabolic co-ordinates, a stream function is obtained for shear flow 
past a finite parabolic cylinder v = v9. The separation constant used to solve | 
this problem takes continuous values from zero to infinity. In order to keep 
the integrals involved convergent, the parabolic cylinder is restricted in extent | 
to the range 0 <u<A. The flow pattern obtained is artificial in that the © 
points « = A on the cylinder are singular points being either sources or sinks. — 
As a result the stream function obtained gives an accurate representation of 
shear flow past a parabolic cylinder only in the vicinity of the nose, which can 
be extended by taking A sufficiently large. | 

The main features of the above shear fields are determined, particular 
attention being paid to the location of the stagnation points. A feature of this 
study is the repeated appearance of the non-dimensional number N = U/l w 
where U and / are a typical velocity and length in the field respectively and w — 
is the vorticity. This number also plays an important part in constant rotational 
flow through a channel. If the latter has straight parallel sides, for example, it 
is easily shown that a necessary condition for no back flow in the channel is 
N = 1/2, where U is the undisturbed velocity along the centre line and / is the 
breadth of the channel. In fact, N which also appears in the force coefficients 
seems to have an importance in rotational flow comparable to the REYNOLD’s 
number in viscous flow. 

In conclusion, it is not suggested that elliptic and parabolic co-ordinates 
exhaust the possible cylinder cross sections for which constant shear flow 
stream functions can be obtained. These co-ordinate systems do, however, 
serve to illustrate the methods described in the present paper. 
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Zusammenfassung 


Tsrens Methode, die Stromfunktion fiir Strémung mit konstanter Drehung 


an einem Zylinder vorbei festzustellen, wird erweitert, um Gebrauch von den 


. 


 natiirlichen Koordinaten des Systems zu machen. Das modifizierte Verfahren 


wird gebraucht, um die Stromfunktionen fiir Stromung mit Drehung an ellip- 
tischen und parabolischen Zylindern und an einer ebenen Platte vorbei festzu- 


stellen. Die Hauptziige der Str6mungsfelder werden untersucht; man richtet be- 


‘sondere Aufmerksamkeit auf die Positionen der Staupunkte. Als ein dausserst 


wichtiger Parameter in dieser Studie erweist sich die dimensionslose Zahl 


_N= U/ol, wo U und l respektiv eine typische Geschwindigkeit und Lange im 


Felde sind und w die Drehung ist. 


(Received: Mai 17, 1954.) 
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Uber die asymptotische Darstellung der Aufspaltung von 
Paaren benachbarter Eigenwerte der Differentialgleichung der 
Spharoidfunktionen’) 


Von Kraus KRICKEBERG, Wiirzburg?) 


Es ist wohlbekannt, dass die Eigenwerte Any?) der Differentialgleichung der 
Spharoidfunktionen 


fa — 24) yan’ + [at 2-24) - 4] 7) = 9 (1 


wobein, m= 0,1, 2,... und m™ < n, sowohl fiir grosse positive 7? (das heisst y >1) 
als auch fiir grosse negative y® (das heisst y* = — i y > 1) asymptotische Ent- 
wicklungen nach den Potenzen y?, y, 1, y7}, ... besitzen’). Im ersten Fall (y > 1) 
haben verschiedene Eigenwerte auch verschiedene derartige Entwicklungen, wah- 
rend im zweiten Fall (y* > 1) die Entwicklungen nur von den Zahlen m und ¢ 
abhangen, wobei g= ”, wenn ” — m ungerade, und ¢g = n +1, wenn »—™ 
gerade. Bei gegebenem m und p (mit p= 0, 1, ...) haben also die zum,=m+4+ 2p+1 
und n,= m+ 2p gehérenden Eigenwerte Ag A= pe) und A (— y**), von denen 


iibrigens der erste einer ungeraden Eigenfunktion Psp (23 — y**) und der zweite 


einer geraden Eigenfunktion Psy, (23 — y*2) entspricht, die gleiche asymptotische 


1) The preparation of this note containing partial results, was sponsored by the European 
Office Air Research and Development Command, U.S. Air Force, Project No. AF 61 (514)-443. 


2) Mathematisches Institut der Universitat. ; 
3) Man findet sie in dem kiirzlich erschienenen Buch von J. Merxner und F.W.SCHAFKE,, 


Mathieusche Funktionen und Sphdaroidfunktionen (Berlin, Géttingen und Heidelberg 1954), dessen 
Bezeichnungen auch hier verwendet werden, auf den Seiten 243 und 246. 
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41 = Ant p= pe hiss Mag yes 
gibt es jedoch eine asymptotische Darstellung nach den Funktionen 
Ba yey ne ee 
b! (m+ p)! 


die sich zum Beispiel ftir die numerische Berechnung von Tafeln der Eigenwerte — 
als niitzlich erweist. Diese Darstellung lasst sich nun nach einem von Herrn 
Professor Dr. J. MEIXNER?) angeregten Verfahren um ein Glied, das heisst bis 
zu k = 2, fortfiihren, und zwar ergibt sich fiir grosse positive y*: 


Entwicklung fiir y* > 1. Fiir die « Aufspaltung », das heisst die Differenz 
y*-k mit k = 0,13), | 

: 

| 


2 (4 yi tl o—2y* 1 . 
aig CE [top een hs TP | 
1 9¢ qe ares ial (oles), q(m? — 1) 2) 
oy esr | 8 ots 2 — 34 q 4 A 2 Ss | 
ees 
{ 2 tee 
rm 1+ 8 {7 


Da die zur Formel (2) fiihrenden Rechnungen recht umfangreich sind, soll 
das Verfahren nur kurz beschrieben werden. Es beginnt mit dem aus MEIXNERS 
und ScHAFKEsS Buch) folgenden, zumindest fiir 0 < 6 Sz <1 richtigen Ansatz 


DS (eat hi NN Sean Bip) a, el eet (3) 


in dem N eine nur von m und ” abhangende Konstante und M q2, mig eine Whitta- 
kersche Funktion bedeutet. Aus (1) und der Differentialgleichung der Whittaker- 
schen Funktionen erhalt man die folgende Differentialgleichung fiir die Funktion 
S) 
hey? ER EE te (ER 2 gE 4 me 1) E44 3S" 8 SR a ae ee Be 
£/2 £2 &’2 &2 
= wage a t 
arti eae RT, ins 


= 0, 


Zusammen mit der Bedingung &(1) = 0 (fiir jedes y*) gestattet sie, die Koeffi- 
zienten der Darstellungen 


Am~c_yy* + C+ Cc, p*-!, (4) 
E~ Dy + Got OL y* + py yt? (5) 


vermittels sukzessiver Bildung von expliziten Darstellungen fiir die Ableitungen 
yj, und deren Integration zu berechnen. Die Zahlen c_,, Cy und c, und die Funktio- 
nen @_y, Po und y, waren bereits bekannt*); der Vergleich jener nach einem 


1) A.a.O., S.319, Formel (14). 

*) Auch an dieser Stelle méchte ich Herrn Mr1xNer noch einmal fiir diese Anregungen danken. 
Ferner danke ich Herrn H. Brunke fiir die Kontrolle eines grossen Teils der Rechnungen. 

Avan On Sao lene) 

4) A.a.O., S. 246, Satz 12, und S. 318. In dem auf S. 318 in der Formel fir € enthaltenen 


Ausdruck fiir m,, das heisst im Koeffizienten von y*—1, muss der letzte Term der ersten eckigen 
Klammer nicht —2, sondern —3/2 heissen. 


A aes ia! 
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anderen Verfahren gewonnenen C_,, Co, C¢, mit den wie eben ermittelten bildet eine 
Kontrolle der Rechnungen. Dariiber hinaus ergibt sich 


‘ gq 1 1 q(1 — m? 
0) flap z|+ rd 


a 
il AL i i 1+ 2 3 : 
x|s=7 n q 
: Gee fee a a a gO ac 8 
; . 1 (eg tier us 3 ‘ 3 1+2 
fies CAs eo 2S a iy 2 
; ; 1 ao 5 
oh =— 
Dm) ae (1 2 ral 
Ist nm, — m ungerade, so genugt PSs, der « Eigenwertbedingung » 
PSO, = De Ds (6) 
wahrend bei geradem , — m gilt 
ops” 
a ye 4 RBC 
elo — 72) = 0. (7) 


Diese Bedingungen lassen sich nun formal auch mit Hilfe der Ansatze (3) und (5) 
erfiillen, wenn man q einmal durch q,(y*) = ¢ + 69,(y*), das andere Mal durch 
qo(y*) = 7 + Og2(y*) ersetzt, wobei 6g, und 6gy in geeigneter Weise zu bestim- 

_ mende Funktionen von y* bedeuten. An die Stelle von p sind dabei die auf Grund 

von 2 p= q — m— 1sinngemass gebildeten Ausdriicke p, bzw. p, zu setzen. Man 
verwendet am bequemsten in (3) statt Mojo, mjg auf Grund von 


MM, /2, m/2 (—#) a ner De eel? F(—?, m + 1 —&) 


_ die konfluente hypergeometrische Funktion F = ,F;, operiert also direkt mit der 
Darstellung aus dem Werk MEIxNeERS und ScHAFKes!). Mit der Abktirzung 
£, = &(0) geht dann die Bedingung (6) tiber in 

F(—p,, m + 1, —&) = 0, (8) 

und analog wird aus (7): 

F(—Pe, m + 1, —&,) [€”(0) &(0)? — (m+ 1) So* + 1 

? 4+ 2F'(—po,m-+ 1, —&) = 9, 
wobei F’ die partielle Ableitung von F nach dem letzten Argument bedeutet. 
Dabei ist also g in den aus (5) resultierenden Darstellungen von ép, &’(0) und €”(0) 
im ersten Fall durch q, und im zweiten durch 4, zu ersetzen. Verwendet man die be- 
kannte asymptotische Entwicklung von F(—p,,m-+ 1, —& )/I'(m-+ 1)*) und ersetzt 
den darin vorkommenden Term 1/I'(—f,) durch (—1)+1 I'(p, + 1) a~'sin (x 6p,), 
so lassen sich in dem Ansatz 

2 (4 p*)t e297" 
0g, ™~ 
41 ~ “pl (mm + 2)! 


1) A.a. O., S. 318, (12). 
2) Vgl. zum Beispiel H. Bucuuorz, Die konfluente hypergeometrische Funktion (Berlin, 


Géttingen und Heidelberg 1953), S. 93, (6); man beachte, dass aus y*>1 folgt —&>1. 


(9) 


[1 + a, y*-* + as y*~*] (10) 
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die Koeffizienten a, und a, nacheinander so bestimmen, dass die ersten drei 


Glieder der resultierenden asymptotischen Darstellung von F(— p;, m + 1, — &) 
verschwinden, dass also (8) in den ersten drei Gliedern erfiillt ist. Durch einen ent- 


sprechenden Ansatz fiir dq, kann (9) befriedigt werden, und dabei zeigt sich 
iibrigens, dass — dq, die gleiche Darstellung (10) wie 6g, hat. Berechnet man 
schliesslich die Reihe (4) fiir 4 einmal mit g, und einmal mit q, statt g, so liefert _ 


die Differenz die gesuchte Formel (2) fiir Ad. 
Das Verfahren kénnte, wenn auch mit betrachtlichem Aufwand, in gleicher 
Weise fortgesetzt werden. Die Bestimmung jedes weiteren Gliedes der Darstellung 


von AA erfordert die Berechnung einer weiteren der in (5) auftretenden Funktio-: 


NED Px. 

Die folgende Zahlentafel dient dem Vergleich einiger direkt berechneten 
Werte von AA mit den entsprechenden durch die rechte Seite von (2) gegebenen 
Werten (AA()), Daneben sind auch noch die Werte aufgefiihrt, die die nur bis zum 
Glied mit y*-1 fortgefiihrte rechte Seite von (2) liefert (4A)1). Es ist f= y*~1. 


t AA AA) 
ae 0,5 1,0892 22, 
ae 0,3 0,30062 0,31675 
0,1 0.5921 (— 5) 0,5936 (—5) 
ar 0,10 0,57629 (—2) 0,52765 (— 2) 
Q 3 . 0,08 0,109385 (— 3) 0,104159 (— 3) 
0,065 0,86 (26) 0,84 (— 6) 
, 1 0,250 1,03934 1,03055 
- 2 ) 0,15 0,38777  (—1) 0,38391 (—1) 
0,10 0,19890 (— 3) 0,19787 (— 3) 
yaa 0,075 OS5318, (— 3) 0,45275 (—3) 
0,060 0:254  {—5) 0,229. (—5) 
fa 0,055 0,20 (6) OAS onde 6) 
AA AAA A= AR 
Peter l 1,0623 — 0,0830 + 0,0269 
f= | 0,30911 — 0,01613 — 0,00849 
0,5924 (—5) = 00015 (—5)" 0000s w= 5s) 
pede 0,57976 (— 2) + 0,04864 (—2) —0,00347 (—2) 
ae 0,109645 (— 3) + 0,005226 (—3) — 0,000260 (— 3) 
0,87 (— 6) + 0,02 (2s 00d (— 6) 
pan 1,11642 + 0,00879 — 0,07708 
aoe 0,39082 (—1) + 0,00386 (—1) —0,00305 (—1) 
0,19919 (— 3) + 0,00103 (—3) —0,00029 (—3) 
aan 0,55940 (— 3) + 0,10043 (—3) —0,00622 (—3) 
Tell O255. {=5) + 0,025 (—5) —0,001  (—5) 
0,20 (— 6) + 0,02 (— 6) 0,00 (— 6) 


1) Ich danke Herrn ABrAmowrrz, National Bureau of Standards, Washington, fiir die Mittei- 
lung der Werte von AA und meiner Frau fiir die Berechnung der Werte von 4A@ una AA), 


y 
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Summary 


Let Ad be the difference of a pair of adjacent eigen-values of the differential 
equation (1) for the spheroidal functions. Until now only the two first terms of 
the asymptotic expansion (2) of 4A for large y*? = — y2 had been known. In this 
note the next term, i. e. the coefficient of y*~2, is given, and the way of calcu- 
lating it is described. 


(Eingegangen: 31. Dezember 1954.) 


On Conditions for Stability of Solutions of 
Pendulum-Type Equations’) 


By James C. LiLo and GEORGE SEIFERT, Lincoln, Nebraska, U.S. A.) 


Introduction 


We consider the equation 
a3 ; - a6 
§+ 0/0) 6=e(0) (6= 4), (1) 


where f and g are twice-differentiable functions of period 2 2 in 8, Ue) coe) tieye 
all 6, g(0) has a finite number of simple zeros in any interval of length 2 2, 


2a 
[g(0) d8>0, 
0 


and « is a positive parameter. If we put y = §, then equation (1) is equivalent to 
the system, ; 
p= e(0)—af()y, O=¥ 


from which we obtain the equation of the phase trajectories in the cartesian 


(6, v) phase plane, 
dy _ g(9) — «f(8) 918) : 
do (0) 

It is known that if equation (2) has a solution of period 2 2 in 9, this solution 
is unique and nonnegative (3) Lets also known that the nonexistence or existence 
of such a solution depends on whether or not « exceeds a critical value «,> 0. 
Furthermore, if no such solution exists, each solution 6(t) of equation (1) is such 
that jim, 6(¢) exists and is finite; i.e., all solutions are stable with respect to a 


singular point of equation (2). ‘ 

The purpose of this note is to determine bounds on @. 

For the special case g(6) = B — sin 0, 0< 6 < 1, and f(@) =1, the best bounds 
to date have been given by HAYEs [2] and Boum [3]. Each gives as an upper 


1) This research was supported by the United States Air Force through the Office of Scientific 


Research of the Air Research and Development Command. 
2) Department of Mathematics and Astronomy, University of Nebraska. 


3) Numbers in prackets refer to References, page 229. 
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bound on «, the value 2 sin(6@,/2) where 9, = arcsin 8. BOHM’s lower bound, 


M4 


sin 69 (5 


-1/2 
sin 65 + > cos a) . 


is greater than Haves’ 
[(3 cos? @) + 1)! —'2 cos @,}!# 


for values of 6, near zero, while the opposite is true for 0, near 2/2. 

Conditions on / and g in addition to conditions on « are known under which 
equation (2) will have a positive solution of period 2 2 in 6 [4]. A study of the 
variation of such a periodic solution in terms of variations in « has given a method 
whereby an upper bound on «, can be determined in terms of a known lower 
bound. However, in the special case mentioned above the values obtained are 
inferior to those given by B6um and Haves especially near B = 1. 

Let us denote by .6;, 0;< O41 7= 0, +1, +2,..., the zeros of g(0). These 
zeros are Clearly the singular points of equation (2); the points (6;, 0) for which 
g’(0,;) > 0 are saddle points, the others being either nodes or focal points. The 
fact that 


27 
|g(0) da > 0 
0 
implies the existence of a 0; such that g’(0;) > 0 and that 
6 
[g(s) ds >0 Oo), 
95 


(cf. the lemma in [5]). It is clearly no loss of generality to assume that 6; = 0 for 
some such 7, which we will do henceforth. 


The Lower Bound 
Let us first assume that 6 = 0 is the only zero of g(@) in 0 < @ < 2 2 for which 


6 
[g(s) ds >0 (0 > 6). 
a 


Consider now the solution (6) of equation (2) for which 


piim, (4) =0 and glim, y’(8) > 0. 


Suppose that for some 0 = b, 0> b> 22, we have y(6) = 0. Multiplying equa- 
tion (2) by y(@) and integrating from 6 = 0 to 6 = b, we have clearly 
b b 
| g(8) d0= « [ #8) (0) a0. (3) 


0 1) 


However, from equation (2) we also have that y’(6) y(9) < g(@) for O< 0<b, 
from which we obtain upon integration, 


2 


a 


0 1/ 
y(8)< (2/9 i) (0 0.2 Ba 
0 


a © 


zi : Beat 
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_ Using this inequality in equation (3) we obtain 
b b 6 a2 

[g(0) d0<«f t(6) (2/08 i dé. (4) 
0 0 0 
3 Let us denote by 0; the least zero of g(@) greater than or equal to b. Consider the 
_ Tegion R of the phase plane bounded by the graphs of y = g(@)/«f(@) and y = 0. 

Since the phase trajectories of equation (2) have positive slope in Rk, we must 
have b = 6,_,. Hence 


05 1/2 


j 22 6 2 
[g(0) d0<« [ 7(6) (2/a ds) d@. 
3 0 ) 0 
If we put G = grin, [ €(0) d0, i=1, 2, ..., the left-hand side of this inequality 
can be replaced by G and hence if 
G 
“> = (5) 


Q7 sk 1/2 
19) (2/06 i) d0 
0 0 


we can only conclude that (22) = 0, which implies the existence of a non- 
negative solution of equation (2) of period 2 2 [1]. Hence a» is a lower bound on Bp. 
If 0 = 0 is not the only zero of g(6) in 0 S 6 < 22 for which 


6 
Les) d(s)>0 (06> 68) 
6 


then each such 6 = 0; for which it is true will by an analogous procedure give a 

_ lower bound «; on «,. In general then we take the greatest of the «; as the lower 
bound on «,. 

Because of the nonelementary integral occuring in the denominator of the 
lower bound a, given by (5) above in the special case g(8) = sin 4 + sin (0 — 95) 
where 0 < 0, < 2/2 and /(6) = 1, an overall comparison between a and the 

- lower bounds given by Hayes [2] and Boum [3] is difficult. 

For 6, > 2/2, % approaches a value between (2) and (2/2)1/*; Le., a 18 
inferior to the values of HavEs and Boum. However, if 0) > 0 we observe that 
a((05) > 2/4, while the slopes of BOum’s and Hayes’ bounds approach 2/z and 
1/2 respectively. This indicates that near 6, = 0, the %») given above is greater 

- than the lower bounds of BOHM and HAyEs. 


The Upper Bound 


We define 
ated 1 Re 2 
£1(9) ne sin (0/2) (0 < 6 <4 9) 
while : 
g,(0) = lim g,(6), and g,(2) =, lim g,(9). 


6—0+ O—>22- 


Since g’(0) and g’(2 ) are each finite, these limits exist and are also finite. Hence 
g,(0) is continuous on 0S 0522 and attains its maximum value M in this 


interval. 


ZAMP VI/16 
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Fix A> 0 and consider the graph of y = h(@) = 2 A sin(6/2) in the phase 
plane. If a is such that 
h(0) h’(0) > g(0) — « f(0) A(0) (0<A< 2a), 


then in this interval the phase trajectories of equation (2) crossing y = h(@) 


have slopes less than the slopes of y = h(9) at the points of contact. Clearly, this — 


will be the case whenever 


= M(A)\r 


(9) — 2 A® cos (6/2) 
2 A f(6) 


Suppose that «, > M(A). Since a,< oo (cf. Theorem 2.1 in [5]), the periodic 
solution y,(0) of equation (2) corresponding to «= «, has at least one zero in 
0<6< 22 (cf. Part I in [5]). It follows that y,(0) = y,(2 2) = 0. 

Tf lim y2(0) > h’(0), then there exists a 6 >0 such that 4,(0) > (6) for 


0 < 6< 0. Consider the trajectory of equation (2) containing the point [a, h(a)] 
where 0< a< 0. It is clear that this trajectory must contain the point (0, 0). 
This is impossible since there is only one trajectory emerging from this point 
with positive slope. Hence, we must have 


jim, v¢(@) S40); ie, y(8) << A(8) (0<0< 2a). 


«> max|* 


Now consider the trajectory containing [b, (b)] where 0< b< 2a. It must 
pass through (6,, 0) where 6; is a common zero of y,(8) and g(@). But this is again 
impossible, since the trajectory at (0;, 0) having negative slope is unique. 

Hence, we must conclude that «, = M(A). If M, = min M(A), then a, S My; 
1.e., M, is an upper bound on «,. 

If M, = max g,(6) and m = min f(@), then it is a simple matter to show that 
My S (My/m2)1”. 

In the special case 


g(6) = sin®, + sin(@— 6,), (6) = 1 


the fact that M, is identical with the bound given by Haves and BOum is due 
to the fact that its derivation follows closely the lines of BGum’s approach [3]. 

In either case, however, the bound can ostensibly be improved by the use of 
a more general function /(#) containing several rather than a single parameter, 
yet having the same general properties as the function 2 A sin (6/2). Needless to 
Say, the complications in the corresponding derivation of M, become significant. 
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Zusammenfassung 


Es gibt zu der Gleichung _ 
B+ a f(8) = gf), 
2% 


wo f(8) = f(@ + 22), g(0) = g(@+ 22), f(6) > 0, [ g(®) d®>0 und ¢(6) ein- 


0 
fache Nullstellung hat, einen Wert « = «, mit der Eigenschaft: ist a > «,, so gibt 
es keine Lésung @(¢) dieser Gleichung, fiir die ® = y(®) = y(8+ 2 x) fiir jedes 0 
ist; ist aber « < «,, so gibt es eine solche Lésung, die mitunter nicht negativ ist. 
Im ersten Fall, « >a,. gilt auch fiir jede Lésung &(¢), dass jim B(t) = 8; < oo. 
a . . aaa 
In diesem Artikel werden fiir «, Grenzen gegeben und mit denjenigen verglichen, 
’ die Hayes [2] und Béum [3] im Falle 
g(0) = sind, + sin(@— @,), f(8) =1 (0 <8, < 5) 
erhalten haben. 


(Received: October 18, 1954.) 


Methode zur Bestimmung von Dimpfungsgréssen in Transonic- und 
Uberschallkanilen 


Von Fritz LuiscHEeR, Emmen’) 


1. Einleitung 


Die im allgemeinen tiblichen Methoden zur Bestimmung des Dampfungsmasses 
an einem Windkanalmodell sind: 
a) Ausschwingversuch. — Dabei kann die vorhandene Dampfung aus dem loga- 
rithmischen Dekrement und der Frequenz der auftretenden Eigenschwingung 


ermittelt werden. 
b) Erzwungene Schwingung. — Ermittlung der Dampfung aus Grosse und Pha- 


senlage der erregenden Kraft. 

Diese beiden Methoden lassen sich in Unterschallkanalen sehr gut anwenden. 
In Hochgeschwindigkeits-, Transonic- und Uberschallkanalen treten jedoch sol- 
chen Versuchen gréssere Schwierigkeiten entgegen. Diese sind einerseits bedingt 
durch die grosseren Geschwindigkeiten, anderseits durch die kleineren Abmes- 
sungen der Modelle und deren Baustoff. 

Eine neue Messmethode, die sogenannte « Rotor-Methode», bei der obener- 
wahnte Nachteile eliminiert sind, wurde ausprobiert und mit Erfolg im Transonic- 


Kanal angewendet. 


2. Bezeichnungen 


Massgebend fiir die instationaren Luftkrafte eines Schwingungsvorganges ist 
die als reduzierte Geschwindigkeit bezeichnete Grosse 


V=v/(@/), 


1) Aerodynamische Abteilung des Eidgendssischen Flugzeugwerkes Emmen. 
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worin 


v = Anblasegeschwindigkeit (m/s) 
@ — Kreisirequenz (s—+) 
1 = Bezugsfliigeltiefe (m) 

Die weitern, im folgenden verwendeten Gréssen sind: 
&® = Massentragheitsmoment (mkg/s?) 
Ms = Moment der Tragheitskrafte (mkg) 
My = Moment der Dampfungskrafte (mkg) 
cg = Dampfungsmass = M¢/(ov F b? o) (—) 
o = Luftdichte (kg/s?/m4) 
F = Fligelflache (m?) 
b = Spannweite : (m) 
gy = Amplitude (—) 
k=) Energie (mkg) 
Ea ikeisnun gs (mkg/s) 


3. Vergleich der Verhaltnisse fiir Messungen bei Unterschall- und 
Transonicgeschwindigkeiten 


Die sich aus den Kanalbedingungen ergebenden Zahlenwerte sind in der nach- 
folgenden Tabelle zusammengestellt: 


Modell : I II 

Kanal Unterschall Transonic 
ModellmaBstab 1S 3:100 
Material Holz Leichtmetall 
1 m 0,555 0,15 
v ms i. 50 i 300 
OPT — 1 HATS 
co = 0,1 0,1 
V —~ it 10 1 10 
7) s—t 90 9 2000 200 
Ms/Mg 35 3,5 280 28 


petit 2 pees Be sl) ieee ote Iles oh a il 


Es zeigt sich also, dass im Transonic-Kanal das Verhaltnis Tragheitskrafte : 
Dampfungskrafte um eine ganze Groéssenordnung ungiinstiger ist als im Unter- 
schallkanal. Die grésseren Tragheitskrafte ergeben auch grossere, die Messung be- 
einflussende Sekundireffekte, wie etwa die Lagerreibung, so dass die zu erfassende 
Dampfung nur noch einen kleinen Teil der die Messung beeinflussenden Gréssen 
darstellt und die Genauigkeit der Ergebnisse entsprechend vermindert wird. 

Beim Unterschallmodell liegen die Kreisfrequenzen w in einem versuchstech- 
nisch giinstigen Gebiet, und auch die Federkrafte beim Ausschwingversuch bzw. 
Erregungskrafte bei der Methode b) lassen sich ohne weiteres aufbringen. Beim 
Transonic-Modell mit den hohen Werten von @ treten auch von dieser Seite her 
Schwierigkeiten in Erscheinung. Dies ist insbesondere bei der Ausschwingmethode 
der Fall, wo fiir eine bestimmte Machsche Zahl die Variation der reduzierten Ge- 


a? tae 
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- schwindigkeit V = v/(w/) nur durch Anderung der Frequenz, das heisst des 


"7 


Federsatzes, moglich ist. 
Um diesen Schwierigkeiten zu begegnen, wurde nach einer Messmethode ge- 


sucht, die einfache und prazise Messung sowie zuverlassige Auswertung erlaubt. 


Zur genauen Messung ohne Stérung des Bewegungsablaufes eignen sich zwei- 
fellos die Zeit ¢ und die Frequenz w am besten. Wenn nun zudem noch die Messung 
mit konstanter Amplitude durchgefiihrt werden kann und sich tiber einen ganzen 
Frequenzbereich erstreckt, kénnen die Messungen rasch, prazis und mit einfachen 
Mitteln erfolgen. 

Die vom eidgendssischen Flugzeugwerk Emmen im Transonic-Kanal ange- 


- wandte «Rotormethode» weist diese Merkmale auf. 


4. Rotormethode 


Auch hier wird ein Schwingungsvorgang zugrundegelegt, die Bewegung des 
-Modelles wird jedoch durch einen Rotor mittels einer Exzentertibertragung er- 
zeugt. 

Die Rotationsenergie des Rotors dient zur Aufrechterhaltung der Schwingungs- 
bewegung. Die Drehzahl des Rotors und damit die Frequenz des Modells nehmen 
entsprechend der aufzuwendenden Dampfungs- und Reibungsarbeit stetig ab. 

Die nachfolgende schematische Skizze zeigt die Einrichtung fiir die Bestim- 
mung des Dampfungsmasses um die Querachse eines Halbmodells. 


| Drehachse 
F - Ca a 


Figur 1 
Prinzipskizze zur Rotormethode. 


Mit einer solchen Einrichtung ergibt sich eine Nickbewegung des Modells, die 
eine konstante Amplitude und in erster Naherung einen harmonischen Verlauf des 
Nickwinkels g aufweist. 

Die vom Rotor abgegebene Energie zerfallt in die beiden Anteile 

Leistung durch Reibungskrafte, 
Leistung der dampfenden Luftkrafte, wovon die letztere zur Ermittlung 
des Dampfungsmasses dient. 
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Die Messung beschrankt sich bei bekannten #y (Modell) und #p (Rotor) auf 
die Bestimmung des Verlaufes ¢= w= f(t), womit die Energie Ep des Rotors 
zu jeder Zeit ¢ bekannt ist. ae 
ER = Or oF 


Die Leistung Lp des Rotors fiir die aufzuwendende Reibungs- und Dampfungs- © 


arbeit wird dann: . 
Lr = Or wm? (<) . 


Der Reibungsanteil des ganzen Systems kann aus einem Standversuch v = 0 
ermittelt werden, wobei das Modell durch eine Masse mit aquivalentem Trag- 
heitsmoment #y ersetzt wird, auf welche keine dampfenden Luftkrafte wirken. 
Diese Leistung wird als é 


7) 
[oe Ce (eee 
; bal ( o )o 
bezeichnet. 
Aus dem Leistungsgleichgewicht 


Lael gail 


kann der Leistungsanteil der ddmpfenden Luftkrafte angesetzt werden zu 


to =~" [(2)- (2) 
@ o /o 
woraus sich das Dampfungsmass ohne weiteres bestimmen lasst. 
Eine mit der «Rotor-Methode» gewonnene Kurve 
cg = f(Ma) 


fiir konstante reduzierte Geschwindigkeit V ist aus der nachstehenden graphischen 
Darstellung ersichtlich. 


C. 
gL 
C. 
¢= (Me) V=konstent 
Ma 
tt] er? 
05 08 07 98 G9 10 


Figur 2 
Qualitativer Verlauf der Dampfungskonstanten in Abhingigkeit der Machzahl. 
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Der messtechnisch unerwiinschte Reibungsanteil sollte auf ein Minimum be- 
schrankt werden, um auch in-der Auswertung des Restanteiles Dampfung gent- 
gender Genauigkeit erreichen zu kénnen. Nun ist aber besonders bei hohen Fre- 
_ quenzen ein grosser Teil der mit «Reibung» bezeichneten Leistung auf die Venti- 
_ lationswirkung des Rotors zuriickzufiihren. Durch Unterbringung des Rades in 
einem evakuierten Gehause lasst sich daher ein gewisser Anteil an Reibungs- 
arbeit eliminieren und damit das Verhaltnis 


Reibungsanteil 
Dampfungsanteil 


giinstiger gestalten. 


5. Zusammenfassung 


Um den messtechnischen Schwierigkeiten, welche sich der Anwendung der 
iiblichen Methoden zur Bestimmung des Dampfungsmasses an einem Windkanal- 
modell im Bereiche grosser Anblasegeschwindigkeiten entgegenstellen, zu begeg- 
nen, ist im Eidgendssischen Flugzeugwerk Emmen bei Messungen im Transonic- 
Kanal die in diesem Bericht beschriebene « Rotor-Methode» angewendet worden. 

Diese Versuchsart erlaubt es, das Dampfungsmass ohne jegliche Stérung des 
Bewegungsablaufes durch einfache Messung der Schwingungsfrequenz in Abhan- 
gigkeit der Zeit auch im Bereich der relativ hohen Kreisfrequenzen zuverlassig zu 
bestimmen. 


Summary 


Test difficulties obviate the application of the methods for the determination 
of damping coefficients of wind tunnel models normally used, namely dying-out 
tests or tests with constrained oscillation respectively. 

To prevent these difficulties the ‘“‘rotor method”’ has been developed and suc- 
cessfully tested in a transonic wind tunnel. 

With this method the model carries out an oscillation with constant amplitude, 
which consumes the energy of a rotor wheel. The damping quantity can be deter- 
mined by measuring the oscillation frequency as a function of time without 
influence of the oscillation. By running a test without model but with an equiv- 
alent mass, damping and frictional components can be separated. 


(Eingegangen: 29. September 1954.) 


Some Notes on the Relation Between the Elastic Moduli of Anisotropic 
Cristallites and Quasi-isotropic Polycrystalline Aggregates 


By Hans H. STADELMAIER, Raleigh, N. C., U.S.A.) 


Introduction 


In any attempt to find for a polycrystalline material the mean value of a 
property that is known for its anisotropic constituent, the crystallite, there will 
be included some simplifications. Such simplifications have been applied with 
varying success to the problem of averaging the modulus of elasticity. In Table I 


1) Department of Engineering Research, North Carolina State College. 
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are compiled, by way of review, the mathematical schemes applied to the problem, | 


together with the simplifying assumptions as stated by the authors. 
Table I}) 


Agreement with 


DESEO GOUs observed modulus 


Author Averaging procedure 


ill 1 equal stresses too low but 
ie in all grains better than 3 


Reuss!) 


HUBER and 


= 1 
= E = ——| E(s.,,, 2). dQ f none best 
Scumrp?) 2) (Sex @) ; 
Vv Ee il P Q) dQ equal strains ante 
Si Srvc nats Q5 “in(eee @) in all grains . 


1) Reuss, A., Z. angew. Math. Mech. 9, 49 (1929). 

2) Huser, A., and Scumip, E., Helv. phys. Acta, 7, 620 (1934); Boas, W., and Scumrp, E. 
Helv. phys. Acta 7, 628 (1934). 

3) Voict, W., Lehrbuch der Kristallphysik (Leipzig and Berlin, 1928), Appendix IT. 


The basic assumptions, which have already been discussed by HuBER and 
SCHMID, are, that the crystals be small compared with the dimensions of the 
specimen, that there exist no voids between the crystallites, and that the cohesion 
of the crystallites be upheld by large deformations in a very small transition 
region near the grain boundaries. The need for the latter condition comes from the 
influence of neighbors upon the crystallites. If, in a polycrystalline aggregate, the 
crystallites could deform freely under simple tension, they would do so by varying 
amounts because of the variations in the orientation of the anisotropic crystal- 
lites. Actually, they also exert upon each other, forces other than those which 
transmit the uniaxial stress in the specimen. Thus the stresses may vary in 
magnitude and direction from grain to grain but the assumption of the small 
transition region will furnish crystallites in which the deformation is essentially 
homogeneous. The inhomogeneities are restricted to a region of negligible volume 
and thus do not contribute to the elastic energy. This simplified description, 
inaccurate as it may be, lends itself readily to the problem at hand, and the 
following analysis will show how it is in error?). 


Elastic Energy 


If W* denotes the elastic energy of a crystallite, the elastic energy of the 
whole aggregate follows by summation: 


N 
w* — PULLAR 


j=l 


1) The Sz and c;;, are the elastic parameters in Vorct’s notation. Q is symbolic for the orien- 
tation function (for the c,; and s,;, Q is of the form a? a3 + «3 a2 + «§ at and dQ is the element 
of solid angle so that 2) = 42). 

*) The uniformity of direction in the grain interior and the abrupt discontinuity at the grain 
boundaries, of the ferromagnetic domain patterns resulting from application of uniaxial stress 
[L. J. Dijkstra, and U. M. Martius, Rev. Mod. Phys. 25, 146 (1953)], furnish a striking illustration 
of the validity of the assumed stress distribution. 
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The elastic energy per unit volume then becomes 


* 
ee cm 1 yw, 

Se en 

d 


where W; is the elastic energy per unit volume of the j-th crystallite. The more 
uniform the grain size is, the better the approximation will be. Therefore the 

' condition of uniform size of the crystallites should be added to the assumptions 
stated in the introduction. In cubic metals the elastic energy of the crystallites 
under uniaxial tension is 


2 


hehe vy 
= S34; (S;_ 24) pF 


u 2 
1 , 2 1 , 2 1 U 2 
mes 17 (C;_» 2.) #25 + Q °22; (C54, 2g) Vz + Q °33i (C;4, 2,) 22; 


Sy C49; (Cip» yap) %27 Vyg + C535 (Cp Qu py) Vyj 225 + C515 (Cpa Qaya) Zap Mag 
where 
Qo. = af 03 + 03 of + of of, Qe = BY E+ BS PE + BRB Qy = vive + YE Vet VEY 
Qyap = Qy— Qe— 2B, 
and Qugy and Qgyo follow by cyclic variation of Qu, Qg, and Qy. The «;, B;, y; are 


the direction cosines between the principal axes of deformation and the crystal 
axes. The prime indicates that the elastic parameters are functions of orientation. 


Discussion of Averaging Procedures 


Case 1 and 2. If all stresses X,,; are identical in magnitude and equal to X,, 
the elastic energy of the aggregate becomes 


‘ll 1 1 
Wa a Wi = 7 Xe & Sj Siw On) » 
7 I 
where 


1 , — 1 / 
N 2, ig = o. = a /% (Sjp, 24) d22 . 
j 


_ The summation can be replaced by the integration over all directions because 
the number of crystallites, N, is large. By comparison with the expression for the 
elastic energy of an isotropic medium, 1/2 X3/E, it follows that the average 
modulus of elasticity is E, = 1/5. This result, which corresponds to case 1 of 
Table I, is obtained by using the assumption of REuss, that the stresses are all 
equal. Actually, the stresses will not be the same in all crystallites. Instead, they 
will form a distribution about a mean value. The result derived above is also 
obtained if 1/2 of the mean square of this stress distribution, which is equal to 
the observable uniaxial stress (within the limitations of the assumptions stated), is 
multiplied by 1/N times the sum over the Sixg- Lhus 


1 1 ths ; oer eam 
2W= ad S147 Sin 2,) X35 © 2 Xai N 2 Sj %e5u- (1) 
i] 
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eelioe dite 


Herein the emphasis is shifted from the requirement of equal stresses, which is no t 


longer upheld, to an approximation for the original summation that cannot be 
evaluated otherwise. From the same point of view, case 2 presents itself as a 
different approximation to the elastic energy’), 


1 2 oe 
Lights ae we Xi Ay 
Wi ar Sid "94 = ; oe (2) 
j mes E; (Sip Q,) 
7 


Approximations (1) and (2) can both be viewed as possible methods of replacing 
every Sj,; by an average value which is subsequently factored out. They differ 
only in the way that the average 5,; is formed. To decide which case leads to 
the better approximation for the elastic energy 2 W, the following reasoning 
may be applied. Approximation (1) but not (2) can also be described by replacing 
every X2; by an average X2. Thus the burden of the approximation can be 
assigned to the stresses in (1) but must be carried by the elastic constants in (2). 
From this it will be appreciated that (1) is the superior approximation because 
the relative deviation of X2; from X2is much smaller than the relative deviation 
of the E; from E, when the crystallites are oriented at random. The difference 
between the approximations (1) and (2) is determined by the difference between 
5, and 1/E. Both have been calculated by integration over all directions, and 
values of E derived therefrom are compiled for several cubic metals in reference 2. 
There one finds that 5;; is systematically greater than 1/E and that 1/E is in close 
agreement with 1/E of the bulk metal. The fact that the “‘best’’ approximation 
under the assumptions made in the introduction is greater than 1/E, means that 
the hypothetical conglomeration of semi-free crystallites forms an aggregate 
which is elastically softer than the actual metal, a logical result. This confirms 
in a different way that (1) is the best approximation. The difference between 57, 
[approximation (1)] and 1/E [approximation (2)] results from the skewness of the 
distribution of the values of the sj,; over the crystallographic directions in the 
crystallite (certain directions occur more frequently than others). Approximation 
(2) thus contains a systematic error in the right direction to compensate for the 
discrepancy (resulting from the oversimplified basic assumptions) between 1/E 
and the best theoretical value. The fact that this inherent correction also has the 
right magnitude, is fortuitous, and any implication that the good agreement with 
the observed moduli proves a posteriori that procedure 2 of Table I deserves 
preference over 1, is unjustified. 
Case 3. The strain energy is given by 


al , 2 , 9 , 9 , 
2W= N a (C11; a a. £905 Vyj iz C33; 3; + 2 C195 ¥ 4; Vag 


, : 4 
+ 2 Cons Vg 2ej + 2 C315 25 ¥ 45) : 
The elastic modulus that would make this strain energy expression equal to 
E 2, =2 , 32 liad = = a 
Die [Gt I + HD + ae et Het BY 


1) It will be noted that, to hold with the same rigor as case 1, the expression leading to case 2 
requires equal elastic moduli in all crystallites. This condition, fulfilled only for a fiber texture, was 
overlooked by Huser and Scumrp, because they did not base their method on any general physical 
principle, such as the law of conservation of energy, which is implicit in this discussion. 


i, 
a4 
Fi 
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where yu is Potsson’s ratio, should be smaller than the observed modulus in the 


tre aggregate because none of the basic assumptions have been 
altered. 
When, in formal analogy to case 1, the strains are assumed to be equal in all 


crystallites, one obtains 


ot sud 
OO Og Py gy ee ay 
j j 


If the average ¢;; and Cy, are formed by integration over all directions, they 


_ become 


and 


From c and c, the modulus of elasticity is derived. If, again, this procedure is 
analyzed for its quality as a mathematical approximation, it amounts to repla- 
cing the c;;; and the c,;,; by average values that are subsequently factored out. 
The skewness of the respective distribution functions leads to a c that is too large 
and a c, that is too small. Both errors contribute to a modulus E, that is too 
large, when they are transmitted through the function 

5 (¢ — ¢) (¢ + 2 ex) 


E, = C+ ty : (3) 


If the relative error, which is of the same nature as in E,, would go proportionally 
into ee there would be the same compensation effect as in E,. Instead, with c 
and c, there are introduced two equal relative errors of opposite signs which, 
because of the nature of function (3), more than double the error in Ee This not 
only explains why cop is larger than the observed modulus but also why it is 
more in error than Fy. 


Conclusion 


The basic assumptions stated in the introduction are an oversimplification. 
They should lead to an average modulus of elasticity that is too small, as in case 1. 
The varying degrees of agreement of the calculated moduli with the observed 
modulus of elasticity, are an effect of the simplified basic assumptions combined 
with the quality of the mathematical approximations implicitly employed. 

The error, that is made when the total elastic energy of a polycrystalline 
aggregate is described as the sum of energy contributions from homogeneously 
deformed crystallites, can be calculated from the difference between the average 
value, obtained by integration over all directions, of the elastic parameter sj,, 
and 1/E, where E is the observed elastic modulus of the polycrystalline aggregate. 
May it be noted that the additivity principle for 1/E, applied by WeEERTS?) to 
composite fiber textures, takes on a natural meaning: An addition of two or 
more 1/E values corresponds to approximation (1). An addition of E values 
instead, will obviously not always lend itself to the fortuitous correction that 


improves E,. 


1) Weerts, J., Metallkunde 24, 101 (1933). 
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Zusammenfassung 


Drei Methoden zur Berechnung des mittleren Elastizitatsmoduls vielkristal- 
liner kubischer Metalle werden formuliert als eine Summierung iiber die Bei- 
trage der Kristallite zur elastischen Energie. Dadurch kommt man zu einer 


zwanglosen Erklarung fiir die Abweichungen der berechneten Mittelwerte vom ~ 


beobachteten Modul des Vielkristalls. 


(Received: September 28, 1954.) 
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Stochastic Processes. By J. L. Doos (John Wiley and Sons, New York, ~ 


1953). 654 pp.; $10.00. 

In diesem umfangreichen Werk findet man die vollstandigste und modernste 
Darstellung der Theorie der stochastischen Prozesse. Unter einem solchen ver- 
steht man irgendeinen zeitlichen Prozess, der mit Hilfe von Wahrscheinlichkeits- 
gesetzen beschrieben werden kann. Als wichtige Beispiele aus der Physik k6nnen 
die Brownsche Bewegung oder Diffusionserscheinungen bezeichnet werden. Die 
Theorie dieser Prozesse und ihre mathematische Darstellung beruht einerseits auf 
der mathematischen Masstheorie und anderseits auf fundamentalen Begriffen und 
Satzen der Wahrscheinlichkeitsrechnung. Der Verfasser stellt die wichtigsten 
Grundlagen fiir eine eigentliche Systematik der stochastischen Prozesse zusam- 
men und behandelt nachher einzeln ihre hauptsdchlichsten Klassen, wie Markow- 
Prozesse mit diskretem oder kontinuierlichem Parameter, stationare Prozesse, 


Prozesse mit unabhangigen oder orthogonalen Spriingen, Theorie der «martin- | 


gales» usw. 

Das Buch ist sowohl fiir die Mathematiker als auch fiir die Physiker von 
grossem Wert. Fiir den Mathematiker gibt es eine schéne Darstellung der Be- 
zichungen zwischen der Masstheorie und der Wahrscheinlichkeitsrechnung, dem 
Physiker liefert das Buch die mathematische Darstellung der stochastischen Vor- 
gange. Der Text ist im allgemeinen klar und erschépfend, kleinere Ungenauig- 
keiten in der Formulierung der Voraussetzungen und ausnahmsweise auch bei 
den Beweisen sowie eine erhebliche Anzahl von Druckfehlern wird der Leser 
selber korrigieren. W. Saxer 


«Scientia Electrica». Zeitschrift (Verlag S. Hirzel, Ziirich), Heft 1, Ok- 
tober 1953; sFr. 4.—. 

Seit dem Sommersemester 1949 werden an der Abteilung fiir Elektrotechnik 
der ETH. in regelmassiger Folge Kolloquien iiber « Moderne Probleme der theore- 
tischen und angewandten Elektrotechnik » durchgefiihrt. Neben Hochschulkreisen 
beteiligen sich auch Fachleute aus der Industrie lebhaft an diesen Vortragen. Die 
«Scientia Electrica», deren erstes Heft im Oktober 1953 im Verlag S. Hirzel, 
Ziirich, herausgekommen ist, setzt sich zur Aufgabe, die Vortrage der erwahnten 
Kolloquien auch weitern Interessenkreisen zuganglich zu machen. 

Das erste Heft enthalt die folgenden Beitrage: M. J. O. Strutt, Tvansistoren 
in Endverstdrkerstufen; G. J. EKKERS, Die Feinwanderung an elektrischen Schalt- 
hontakten; O. E. GERBER, Koronaverluste an Hochspannungsleitungen. 


j 


J 
< 

A 
J 
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Die Hefte erscheinen in zwangloser Folge, wobei jeweils 4 bis 6 Stiick in 


_ einen Band zusammengefasst werden. Es ist vorgesehen, neben den Kolloquiums- 


-vortragen auch andere Aufsadtze aus der Industrie und aus Hochschulkreisen auf- 


_ zunehmen, sofern sie in den Rahmen passen. 


Prof. Dr. M. Strutt, der die Kolloquien leitet, zeichnet als Herausgeber der 


«Scientia Electrica». Das Erscheinen dieser auch drucktechnisch vorziiglich aus- 


5 


-gestatteten Zeitschrift kann nur begriisst werden. E, Baumann 


Projective Geometry and Projective Metrics. By HERBERT BUSEMANN 


and Paut J. Ketry (Academic Press Inc., Publishers, New York, 1953). 332 S., 


97 Abb.; $6.00. 
Diese «Projective Geometry» ist kein Lehrbuch im Stil der klassischen Dar- 


 stellung etwa von TH. Reve; gleich zu Beginn werden Koordinaten eingefihrt. 


An Eleganz steht die analytische Methode der synthetischen nicht nach und hat 


-zudem den Vorteil grésserer Allgemeinheit. Zunachst kommen projektive Ebene, 


Polaritaten, Kegelschnitte und affine Geometrie zur Sprache. Eine eingehendere 
Diskussion des metrischen Raumes leitet iiber zur projektiven Metrik und zur 
Minkowskischen, Euklidischen und Hilbertschen Geometrie. Die grundlegenden 


' Begriffe, wie Entfernung, Bewegung, Orthogonalitat, Flacheninhalt, werden sehr 


sorgfaltig eingefiihrt. Ausfiihrlich werden dann die nichteuklidischen Geometrien 
zundchst in der Ebene und dann im Raum behandelt. 

Das Buch, das zahlreiche, gut ausgewahlte Ubungsaufgaben enthalt, ist fiir 
Studierende in den mittleren Semestern geschrieben und vermittelt eine vorziig- 
liche, moderne Einfiihrung in die projektive Geometrie. E. Roth-Desmeules 


Die Gestalt der elektrischen Freileitung. Von MiLan VipMar (Verlag 
Birkhauser, Basel 1952). 200 S., 49 Fig.; sFr. 19.75. 

Das Buch gliedert sich in sieben Kapitel: Gestaltproblem; Spannweite ; 
Querschnitt durch die Dreiphasenleitung; Durchmesser des Ubertragungsleiters ; 
Zustinde des Ubertragungsleiters; Problem des Leitermetalls; Elektromagne- 
tische Gestalt der Freileitung. 

An den Anfang des Buches stellt der Verfasser die Kelvinsche Regel tiber die 
wirtschaftlichste Bemessung von Leiterquerschnitten, welche besagt, dass ein 


_ stromfiihrender Leiter so zu bemessen ist, dass die Kosten der jahrlichen Strom- 


warmeverluste gleich hoch werden wie die Jahreskosten fiir Amortisation und 
Kapitaldienst. Diese Regel wird vom Verfasser fiir Freileitungen ausgebaut; sie 
zieht sich wie ein roter Faden durch das ganze Buch. 

Das 3. Kapitel erwahnt die Begriffe des Wellenwiderstandes und der nattr- 


- lichen Leistung einer Freileitung sowie den Begriff des mittlern geometrischen 


Abstandes der Leiter bei unsymmetrischen Drehstromleitungen, ohne dabei auf 
die Begriindung einzugehen. Erwahnt wird auch die vom Verfasser als « Essener 
Formel» bezeichnete Gleichung fiir den im Hinblick auf Koronaverluste min- 
destens erforderlichen Leiterdurchmesser: d = U,,/9, wobei d = Leiterdurchmes- 
ser in Millimeter, U, = verkettete Betriebsspannung in Kilovolt bedeuten. 
Diese Formel entspricht einer Leiterfeldstarke von 16,7 kV/cm und soll fur 
Netze mit starrer Nullpunktserdung gelten, wahrend der Autor fiir geléschte 
110-kV-Netze U,/8 empfiehlt (14 mm Durchmesser). 110 kV stellt nach den 
Rechnungen des Autors jene Betriebsspannung dar, oberhalb welcher die natur- 


liche Leistung zu steigender Bedeutung kommt und oberhalb welcher das Cu als 


Leitermetall dem Al weichen muss. Dies trifft bei kritischen Leiterdurchmessern 
von 12...14 mm bereits zu. 


~~. x L 
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Im 5. Kapitel geht der Autor einen neuen Weg, um eine Zustandsgleichung 
fiir den Leiter abzuleiten, welche Leiterelastizitat, Temperatur, Durchhang und 
Seilbelastung verkntipfen. 

Im 6. Kapitel iiber die verschiedenen Leitermetalle vertritt der Autor mit 


Vehemenz die Auffassung, dass beim Vergleich von Leitungen aus Cu und Al © 
nicht auf gleichen elektrischen Widerstand (Querschnittsverhaltnis 1: 1,7), sondern 


auf gleiches Gewicht pro Laufmeter abzustellen sei (mechanische Aquivalenz). Das 
bedingt bei der Al-Leitung 3,3fachen Querschnitt des Cu-Leiters. Bei diesem Ver- 
gleich bleiben die mechanischen Beanspruchungen der Masten (abgesehen vom 
Winddruck und von eventuellen veranderten Eislasten bei anderem Seildurch- 
messer) unverdndert. Der elektrische Widerstand ist beim mechanisch A4quiva- 
lenten Al-Leiter kleiner, damit auch seine Jouleschen Verluste. Infolge dieses 
reichlichen Al-Querschnittes betrachtet der Autor alle kiinstlichen Verstarkungen 
der Al-Leiter, zum Beispiel durch Stahlseile usw., als iiberfliissige Fehlentwick- 
lungen. Seiner Auffassung nach sollen Leiter bis 10 mm Durchmesser aus Cu, 
starkere Leiter mit mechanisch 4quivalentem Querschnitt aus Al ohne jede Stahl- 
seele bestehen. 

Das 7. Kapitel ist auf der Vorstellung vom Wanderwellencharakter der Ener- 
gietibertragung aufgebaut. Jeder stationdre Zustand lasst sich bekanntlich aus 
Hin- und Riickwellen zusammensetzen. Diese Tatsache wird vom Autor in ori- 
gineller Weise zur abgekiirzten Ableitung der Leitungsgleichungen beniitzt. Auch 
wird ein Kreisdiagramm der verlustlosen Leitung abgeleitet, das der Verfasser 
nach KoZELJ bezeichnet. 

Das Buch zeigt einen stark persénlichen Charakter. Die Lektiire ist fiir den 
Fachmann ohne Zweifel sehr anregend, auch dann, wenn sich seine Auffassungen 
vielleicht nicht durchweg mit jenen des Verfassers decken. Unter den eher spar- 
lichen Biichern iiber Leitungen stellt es eine Bereicherung dar. K. Berger 


Ferroelectricity. By E. T. JAyNrEs. (Princeton University Press, Princeton 
1953). 136 pp.; $2.-. 

Dieses als erstes einer Reihe «Investigations in Physics» erschienene Bandchen 
verleugnet seine Entstehung als Ausarbeitung einer Dissertation in keiner Weise. 
Die verschiedenen Teile des behandelten Stoffes werden demnach mit unter- 
schiedlicher Griindlichkeit bearbeitet. Die experimentellen Grundlagen nehmen 
bei fleissiger Zitierung von Originalarbeiten einen relativ kleinen Raum ein. Eine 
hiibsche und als Einfiihrung glanzend geeignete Zusammenstellung bisheriger 
Theorien tiber Ferroelektrizitat lauft in eine etwas breitere Darstellung einer vom 
Autor zusammen mit P. E. WIGNER bisher erst in ihren Grundziigen (Phys. Rev. 
79, 213 [1950]) publizierte Behandlungsweise aus. Im Zusammenhang damit 
wird eine von der Ewaldschen Methode verschiedene und in gewissen Fallen vor- 
teilhaftere Berechnungsweise fiir Lorentz-Faktoren angegeben. Das Literatur- 
verzeichnis ist bis Ende 1951 nachgefiihrt. 

Die Art der Darstellung ist klar und anregend, indem nicht nur auf die Erfolge, 
sondern auch auf die Schwéchen im gegenwartigen Verstandnis der ferroelek- 
trischen Erscheinungen hingewiesen wird. H. Labhari 


Molecular Theory of Fluids. By H. S. Green (North-Holland Publishing 
Company, Amsterdam, 1952). virt + 264 pp., 13 Fig.; fl. 20.- (€2/-/-). 

_ Dieser im Rahmen einer Reihe von Monographien iiber Deformation und 
Fliessen herausgegebene Band stellt einen wohlgelungenen Versuch dar, das bis- 
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_ herige theoretische Wissen iiber die Materie im kondensierten, iliessbaren Zustand, 
also als Glas, Fliissigkeit oder dichter Dampf, zu sammeln. Die im Vergleich zu 
den vielen Originalarbeiten leicht lesbare Darstellung geht aus von der Diskus- 
sion, der diesen Zustand charakterisierenden, r6ntgenographisch bestimmbaren 
Dichteverteilung um einen Atom- oder Molekiilschwerpunkt. Einerseits werden 
' die Beziehungen dieser Funktion zu den zwischenmolekularen Kraften und damit 
auch zum zweiten Virialkoeffizienten abgeleitet. Andererseits ergeben sich aus der 
durch einen Geschwindigkeitsgradienten hervorgerufenen Veranderung dieser 
'Dichteverteilung die Fliesseigenschaften. Interessante Abschnitte sind der Be- 
handlung von Kondensations- und Erstarrungsvorgangen gewidmet. Ein letztes 

_ Kapitel enthalt quantenmechanische Verfeinerungen, die bei tiefen Temperaturen 
wichtig werden. In diesem Zusammenhang wird auch das Problem der Super- 

- fluiditat von He II beriihrt. H. Labhart 


Numerical Analysis. By D. R. HartreE (Oxford University Press, Lon- 

don, 1952). 287 pp., 31 figs.; 30s. 

Das vorliegende Werk gibt einen Uberblick iiber das ganze Gebiet der prak- 
tischen Analysis. Beginnend mit einigen Bemerkungen iiber Rechentechnik im 
_allgemeinen und iiber den Gebrauch von Rechenmaschinen, behandelt es Diffe- 

renzenrechnung, Interpolation und numerische Integration. Weiter folgen lineare 
Gleichungssysteme, algebraische Gleichungen sowie partielle Differentialglei- 
chungen (insbesondere Relaxation). Das Werk schliesst mit einem Kapitel tiber — 
_ programmgesteuerte Rechenmaschinen, in welchem dem Leser gezeigt wird, wie 
man ein Problem mit einer elektronischen Rechenmaschine lost. Es darf hervor- 
gehoben werden, dass sich dieses Lehrbuch nicht mit der tiblichen Behandlung 
der numerischen Methoden begniigt, sondern auch auf die Besonderheiten und 
Nachteile derselben hinweist, was besonders fiir den unerfahrenen Leser sehr 
wertvoll sein diirfte. _ HA. Rutishauser 


Tratado de Nomo¢rafia. Por J. C. BELGRANO (Instituto técnico de la con- 
struccion y del cemento, Madrid 1953), 387 pag.; 125 ptas. 

Dieses umfangreiche Handbuch der Nomographie wird durch ein Vorwort von 
Prof. J. Rey Pastor eingeleitet und erwahnt A. Lopez NIETO und J. M. UrR- 
CELAY als Mitarbeiter. Es vermittelt einen umfassenden Uberblick tiber die in 
der praktischen Nomographie angewandten Methoden. In zwei Hauptkapiteln 
sind die Netztafeln und die Fluchtentafeln fiir Beziehungen zwischen drei und 
mehr Variablen behandelt. Inhaltsreich ist auch der Abschnitt tiber Nomo- 

gramme mit transparenten beweglichen Blattern. Das bemerkenswerte Werk 
enthalt sehr viele Anwendungsbeispiele, die, dem Herkommen des Verfassers- 
entsprechend, vorwiegend technischer Natur sind. Daneben finden auch allgemeine 
theoretische Betrachtungen den gebiihrenden Platz. Zu erwahnen ist die einge- 
hende Beriicksichtigung von Beziehungen mit vielen Veradnderlichen, die in der 
Technik ja haufig genug angetroffen werden. E. Vollm 


Sample Survey Methods and Theory. By M. H. Hansen, W.N. Hur- 
witz, W. G. Mapow (John Wiley and Sons, New York, 1953). Volume 1: Methods 
and Applications, 638 pp.; $8.-. Volume 2: Theory, 332 pp.; $8.—. 

Diese zwei Bande kénnen als eigentliches Handbuch oder als ein Standard- 
werk sowohl in theoretischem als auch in praktischem Sinne fiir die Anwendung 
des Stichprobeverfahrens bei statistischen Aufnahmen bezeichnet werden. Die 
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Verfasser, amerikanische Statistiker, angestellt beim Volkszahlungsbiiro ee 
USA. bzw. in einem statistischen Forschungslaboratorium, stehen bei ihren | 
statistischen Untersuchungen in einem grossen Kontinent stets vor der funda a 
mentalen Frage, wie man durch zweckmassige Wahl einer Stichprobe auf die 
Allgemeinheit schliessen kénne und mit welcher Sicherheit solche auf Stichproben . 
beruhende Antworten gewertet werden diirfen. Die Verfasser geben im ersten — 
Band eine breit angelegte Darstellung der verschiedenen Auswahlmethoden der _ 
Stichproben und der Antworten samt ihrem Sicherheitsgrad, die man mit ihrer 
Hilfe erhalt. Im zweiten Band wird mittels im wesentlichen elementarer Methoden 
der Wahrscheinlichkeitsrechnung die theoretische Begriindung fiir die im ersten _ 
Band dargestellten Methoden und statistischen Tests gegeben. Zahlreiche Bei- 
‘spiele, Figuren und Aufgaben erlautern den Text. Die beiden Bande diirften vor _ 
allem fiir praktisch arbeitende Statistiker und auch Studierende der Statistik — 


bestimmt sein. Sie kénnen aber fiir die Stichprobentheorie wohl als das voll- 


standigste Werk der Gegenwart betrachtet werden. W. Saxer — | 


Methods of Mathematical Physics, volume I. By R.Courant and © 
D. Hirsert (Interscience Publishers, New York and London, 1953). 561 pp., 
27 figs.; $9.50. 

Dieses allgemein bekannte Standardwerk der mathematischen Methoden der 
Physik, der « Courant-Hilbert», ist nun in einer ersten englischen Auflage er- 
schienen. Ihr liegt die zweite deutsche Auflage zugrunde, doch sind eine Reihe 
von Zusaétzen und Anderungen angebracht worden. Zu einer weitergehenden 
Modernisierung und Umarbeitung reichte allerdings die Zeit nicht. Der Inhalt 
gliedert sich in die folgenden Kapitel: I. The Algebra of Linear Transformations 
and Quadratic Forms; II. Series Expansions of Arbitrary Functions; III. Linear 
Integral Equations; IV. The Calculus of Variations; V. Vibration and Eigenvalue 
' Problems; VI. Application of the Calculus of Variations to Eigenvalue Problems; 
VII. Special Functions Defined by Eigenvalue Problems; Appendix: Transfor- 
mation of Spherical Harmonics (von W. Macnus); Additional Bibliography 
(beriicksichtigt auch die neueste. Literatur). —- Dieses in vorziiglicher Ausstattung 
erschienene Werk wird auch in Zukunft den Mathematikern und Physikern als 
unentbehrlicher Helfer die besten Dienste leisten. E. Roth-Desmeules 


Praktische Mathematik. Von H. von SANDEN (Verlagsgesellschaft B. G. 
Teubner, Leipzig und Stuttgart 1953). 128 S., 25 Abb.; DM 3.20. 

Der bekannte Verfasser legt hier in dritter, erweiterter Auflage sein Bandchen 
Praktische Mathematik vor. Es ist eine vorwiegend fiir Studierende der Ingenieur- 
wissenschaften geschriebene Einfiihrung, diirfte aber auch dem in der Praxis 
stehenden Ingenieur und Techniker gelegentlich von Nutzen sein. Die wichtigsten 
zur Behandlung kommenden Themen sind: Rechenschieber, Funktionspapiere; 
graphische und numerische Integration; Interpolation in einer Tabelle; Lésung 
von Gleichungssystemen nach Newron und durch Iteration; Approximation von 
Funktionen durch Taylor-Reihen, Polynome und trigonometrische Funktionen; 
Glatten von Funktionen; Ausgleichung nach der Fehlerquadratmethode und mit 
Gewichten; Korrelation. —- Die verschiedenen Uberlegungen sind durch gute Bei- . 
spiele aus der Praxis erlautert; ein besonderes Gewicht wird auf Fehlerabschiat- 
zungen gelegt. - Das Buch erreicht das gesteckte Ziel vollauf. 

E. Roth-Desmeules 


ELEKTRONENMIKROSKOPE 
MIKROTOME 
ELEKTRONENDIFFRAKTOGRAPHEN 


HOCHSPANNUNGSOSZILLOGRAPHEN 


TRUB, TAUBER: ZURICH 


Neuerscheinung aus dem Birkhiuser Verlag 


WILLIAM PRAGER 


PROBLEME 
DER PLASTIZITATSTHEORIE 


Von Professor Dr. WILLIAM PRAGER, 
Brown University, Providence, R. I., USA. 


100 Seiten mit 52 Figuren, Ganzleinen Fr. 12.50 (1955) 
(Sammlung «Lebr- und Handbiicher der Ingenieurwissenschaften», Band 17) 


INHALT 


I. Mechanisches Verhalten plastischer Stoffe. II. Mechanisches Verhalten 
plastischer Tragwerke. III. Tragleitverfahren. IV. Endliche plastische 
Formanderungen. 


Das Werk stellt im wesentlichen eine Niederschrift von Gastvorlesungen 
dar, die der Verfasser im November und Dezember 1954 an der Eidgenos- 
sischen Technischen Hochschule in Ziirich hielt. Die Beschrankung auf 
ausgewahlte Probleme der Theorie ideal-plastischer Stoffe erméglicht es 
dem Verfasser, den Leser bis an die Front gegenwartiger Forschung zu fuhren, 
ohne vorherige Bekanntschaft mit der Plastizitatstheorie vorauszusetzen. 
Die Mehrzahl dieser Probleme entstammen der Fachliteratur der letzten 
Jahre und sind im vorliegenden Band wohl zum ersten Male in Buchform 
pehandelt. Die Darstellung zielt hauptsachlich darauf ab, im Leser ein 
intuitives Verstandnis fiir die Besonderheiten des plastischen Verhaltens 
zu erwecken. Zahlreiche Literaturverweise geben dem Leser die Méglich- 
keit, dieses Verstandnis durch weiteres Studium zu vertiefen. 


Zu beziehen durch die Buchhandlungen 


BIRKHAUSER VERLAG - BASEL UND STUTTGART 


ist am 10. Juni 1955 
70 Jahre alt geworden. 


Folgende Werke von ihm sind 
in unserem Verlag erschienen: 


Die 
mathematische Denkweise 


_Dritte Auflage 1952. 122 Seiten mit | 


11 Figuren und Notenbeispielen und 

9 Tafeln. Ganzleinen Fr.14.55. Samm- 

lung «Wissenschaft und Kulturs, 
Band 1. 


Elemente der Philosophie 
und Mathematik 


Eine Anleitung zum inhaltlichen Den- 

ken. 116 Seiten. Ganzleinen Fr.11.45 

(1952). Sammlung «Wissenschaft und 
Kultur», Band 6. 


—= 


Noch 1955 erscheinen: 
Die geistige Arbeit 


Sammelband der von Professor Dr. 
ANDREAS SPEISER in den letzten Jah- 
ren gehaltenen Vortrage, 


Die Theorie der Gruppen 
von endlicher Ordnung 


Vierte erweiterte und erganzte Nach- 
druckauflage des bekannten Werkes. 


————— eee 


BIRKHAUSER VERLAG 
BASEL UND STUTTGART 


| Prof, Dr. Andreas § "peiser 


Soeben erschienen! — 


Einfihrung in die 


Quantenelektrodynamik — 


Von asin 
Dr. WALTER E. THIrRING . 
Institut fiir theoretische Physik . 
an der Universitat Bern 
1955. XII und 122 Seiten, 18 Di 
Ganzleinen. DM/sFr. 17.50 


Das Buch will den gesichertsten Teil der 
quantisierten Feldtheorie, namlich die Quan- 
tenelektrodynamik, in ihren wesentlichsten 


' Punkten darstellen. Dabei wurde getrachtet, 


méglichst alles zu bringen, was zum physika- 
lischen Verstandnis nétig erscheint, und lie- 


ber auf manche mathematisch-formalen Ein- — 


zelheiten zu verzichten. Das Buch soll kein 
Handbuch, sondern eher eine Sammlung der 
gesicherten, befriedigenden und interessanten 
Resultate sein, die man bisher gewonnen hat. 
Fiir die Lektiire des Werkes (das aus Vor- 
lesungen entstand, die der Autor am Institute 
of Advanced Studies in Dublin gehalten hat), 
ist Analysis und lineare Algebra erforderlich, 
weniger bekanntes mathematisches Riistzeug 
ist im Anhang behandelt; an ph i 
Vorkenntnissen ist spezielle Relativitats- 
theorie und elementare Quantenmechanik 
vorausgesetzt. * 
Ein Buch, das allen Professoren und Do- 
zenten eine wertvolle Grundlage fiir ihre Vor- 


lesungen gibt, und dariiber hinaus es Physi- - 


kern und fortgeschrittenen Studenten ermdg- 
licht, sich in die moderne Entwicklung einzu- 
arbeiten und mit der modernen Grundlagen- 
forschung vertraut zu machen, 7 


Die Laplace-Transformation 
und ihre Anwendung 


Von : 
P, Funx, H. Sacan und F.Szuie - 


Professor, Dozent an der Technischen 
Hochschule Wien 


1953. VII und ro6 Seiten, 18 Diagramme, 
kartoniert DM/sFr. ro.— 


Das Buch schlagt die Briicke von der rein 
mathematischen Literatur zum Techniker, 
dem mit diesem Hilfsmittel ein wirkungsvol- 
les Instrument zur Lésung praktischer Pro- 
bleme in die Hand gegeben ist: «Fiir einen 
nicht an rein mathematischen, sondern prak- 
tischen Problemen Interessierten stellt diese 
Verdffentlichung den bequemsten Zugang in 


deutscher Sprache dar,» (Physikalische Blat- 


ter, 1954, Nr. 5). 

s.., Diese vermittelnde Darstellung ist vor- 
zuiglich gelungen und sehr zu begriissen...» 
(ZAMP, 1955, Nr. 1). : 


Einzelprospekte stehen zur Verfiigung 
In jeder Buchhandlung 


VERLAG FRANZ DEUTICKE, WIEN 


a | 


